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ABSTRACT 

We discuss the Galactic foreground emission between 20 and 100 GHz based on observations by Planck and WMAP. The total intensity in this 
part of the spectrum is dominated by free-free and spinning dust emission, whereas the polarized intensity is dominated by synchrotron emission. 
The Commander component-separation tool has been used to separate the various astrophysieal proeesses in total intensity. Comparison with radio 
recombination line templates verifies the recovery of the free-free emission along the Galactic plane. Comparison of the high-latitude Her emission 
with our free-free map shows residuals that correlate with dust optical depth, consistent with a fraction 30 %) of Her having been scattered by 
high-latitude dust. We highlight a number of diffuse spinning dust morphologieal features at high latitude. There is substantial spatial variation 
in the spinning dust speetrum, with the emission peak (in f) ranging from below 20 GHz to more than 50 GHz. There is a strong tendeney for 
the spinning dust component near many prominent H ii regions to have a higher peak frequency, suggesting that this increase in peak frequency is 
associated with dust in the photo-dissociation regions around the nebulae. The emissivity of spinning dust in these diffuse regions is of the same 
order as previous deteetions in the literature. Over the entire sky, the Commander solution finds more anomalous microwave emission (AME) than 
the WMAP component maps, at the expense of synchrotron and free-free emission. This can be explained by the difficulty in separating multiple 
broadband components with a limited number of frequency maps. Future surveys, particularly at 5-20 GHz, will greatly improve the separation by 
constraining the synchrotron spectrum. We eombine Planck and WMAP data to make the highest signal-to-noise ratio maps yet of the intensity of 
the all-sky polarized synchrotron emission at frequencies above a few GHz. Most of the high-latitude polarized emission is associated with distinct 
large-scale loops and spurs, and we re-discuss their structure. We argue that nearly all the emission at 40° > / > -90° is part of the Loop I strueture, 
and show that the emission extends much further in to the southern Galactic hemisphere than previously recognised, giving Loop I an ovoid rather 
than cireular outline. However, it does not eontinue as far as the ''Fermi bubble/microwave haze”, making it less probable that these are part of 
the same structure. We identify a number of new faint features in the polarized sky, including a dearth of polarized synchrotron emission directly 
correlated with a narrow, roughly 20° long filament seen in Her at high Galaetie latitude. Finally, we look for evidence of polarized AME, however 
many AME regions are significantly contaminated by polarized synchrotron emission, and we find a 2cr upper limit of 1.6 % in the Perseus region. 

Key words. Diffuse radiation - ISM: general - Radiation mechanisms: general - Radio continuum: ISM - Polarization 
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1. Introduction 

Diffuse Galactic radio emission consists of a number of distinct 
components that emit via different emission mechanisms (no¬ 
tably synchrotron, free-free, and spinning dust). There is con¬ 
siderable interest in understanding these components in order 
to subtract foregrounds cleanly from cosmic microwave back¬ 
ground (CMB) data and as a probe of the physics of the interstel¬ 
lar medium (ISM) and Galactic structure. The separation of the 
diffuse foregrounds into their separate constituent components is 
extremely difficult, since their spectra (except for free-free emis¬ 
sion) are not well-known, and they have comparable intensities 
at microwave frequencies (Leach et al. 2008). 

At low frequencies (< 10 GHz), synchrotron radiation from 
electrons spiralling in the Galactic magnetic field dominates the 
sky. At frequencies around 1 GHz the spectral index (T oc yP) 
P ^ -2.7, while at higher frequencies it appears to steepen 
to p j^-3.0. This broadly agrees with theoretical calculations, 
such as those from the GALPROP^ code (Orlando & Strong 

2013) ; however, a detailed comparison has not been made up 
to now, due to difficulties in the component separation process. 
Synchrotron radiation is intrinsically highly polarized and is a 
strong polarized CMB foreground up to around 100 GHz. Free- 
free emission emits over a range of radio frequencies and can 
be significant up to about 100 GHz owing to its flatter spec¬ 
tral index {p = -2.\). Anomalous microwave emission (AME) 
is an additional foreground that has been detected at frequencies 
10-60 GHz (Leitch et al. 1997; de Oliveira-Costa et al. 2004; 
Finkbeiner 2004; Davies et al. 2006; Planck Collaboration XX 
2011; Planck Collaboration XXI 2011; Planck Collaboration 
Int. XV 2014) above the expected levels of synchrotron and 
free-free emission. The most plausible origin for AME is elec¬ 
tric dipole radiation from tiny, rapidly spinning dust grains 
(Draine & Lazarian 1998), which provides an excellent fit to the 
data, particularly from well-studied molecular clouds (Planck 
Collaboration XX 2011). However, the diffuse emission at high 
latitudes has still to be definitively identified, and other possibil¬ 
ities exist, such as magneto-dipole radiation from thermal fluctu¬ 
ations of magnetized grains (Draine & Lazarian 1999; Liu et al. 

2014) . The polarization of AME is of great interest to CMB cos- 
mologists, since a significant level of polarization could hamper 
measurements of B modes. Current measurements place upper 
limits on the polarization of AME at a few per cent or below 
(Dickinson et al. 2011; Macellari et al. 2011; Rubino-Martm 
et al. 2012a). 

In the 2013 Planck^ results based on spectral fitting, the com¬ 
ponent separation procedure did not achieve a separation of the 
distinct low-frequency components, instead lumping these into 
a single low-frequency foreground (Planck Collaboration XII 
2014). However, other methods have been used to do this under 
specific assumptions; Planck Collaboration Int. XII (2013) used 
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(USA). 


correlated component analysis in the southern Gould Belt re¬ 
gion where synchrotron emission is relatively smooth and weak, 
while Planck Collaboration Int. XXII (2015) used template fit¬ 
ting to separate components that were correlated with each spa¬ 
tial template. In the 2015 analysis (Planck Collaboration X 2016) 
we include the WMAP (Bennett et al. 2013a) and 408 MHz 
(Haslam et al. 1982) all-sky maps to enable a separation of syn¬ 
chrotron, free-free, and AME. The separation is based on fitting 
a spectral model to each sky pixel using Commander, a Gibbs 
sampling code (Eriksen et al. 2008). 

The current Planck data release also includes polarization 
information. At angular scales larger than 1° the polarized emis¬ 
sion is dominated by synchrotron and thermal dust radiation. 
Synchrotron radiation is more important at low frequencies and 
dominates the polarized Planck maps at 30 and 44 GHz. This 
synchrotron emission comes mainly from the Galactic plane and 
also from filamentary structures that can extend over 100° across 
the sky. While these “loops” and “spurs” have long been known 
in total intensity (Quigley & Haslam 1965; Berkhuijsen et al. 
1971; Berkhuijsen 1971), their pattern is clarified by the new 
polarization maps. Vidal et al. (2015) present an analysis of the 
filaments using WMAP data. They have catalogued them, stud¬ 
ied the polarization angle distribution and polarized spectral in¬ 
dices, and presented a model to explain the origin of some of 
them. Here, with the new Planck polarization data, we expand on 
this, identifying new polarized features and understanding better 
the already known ones, thanks to the improved signal-to-noise 
(S/N) ratio of the data. 

We begin in Sect. 2 by summarizing the data sets employed, 
including a brief discussion of some of the most important sys¬ 
tematic effects. We then take a first look at the intensity fore¬ 
grounds in Sect. 3 by employing a constrained internal linear 
combination (ILC) algorithm to remove the CMB and free-free 
emission components whose spectra are well known. In Sect. 4 
we then discuss each of the low frequency foreground compo¬ 
nents (AME, free-free, and synchrotron emission) as determined 
by the Coinmander algorithm. We use these new maps to investi¬ 
gate the distribution and spectra of these foregrounds. In Sect. 5 
we use ^MAPjPlanck polarization maps to study the distribu¬ 
tion of the low frequency polarized foregrounds. We make a high 
S/N ratio synchrotron map by combining the WMAP/P/a^ck 
data into a single product, assuming a power-law model for syn¬ 
chrotron emission. We discuss in detail the large-scale features 
in the low frequency polarized sky including the well-known 
loops and spurs. Other results will be published in separate ar¬ 
ticles, including the diffuse synchrotron power spectrum and 
further analysis of AME polarization constraints from Galactic 
clouds. We conclude in Sect. 6. 

2. Data 

Table 1 fists the primary radio data sets used in our analysis, with 
properties that are important here. We note that brightness tem¬ 
peratures are in the Rayleigh-Jeans convention unless otherwise 
stated (WMAP and Planck data are natively in CMB thermody¬ 
namic temperature). We now briefly discuss these data sets and 
particular issues that are relevant to our analysis. 

2.1. WMAP and Planck data 

The primary data used in the analysis are the 9-year WMAP 
maps (Bennett et al. 2013a) and the 2015 release of Planck 
maps (Planck Collaboration I 2016; Planck Collaboration II 
2016; Planck Collaboration VIII 2016) from the Low Erequency 
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Table 1. Summary of radio data sets used, with selected properties. 


Map 

Instrument 

^0 

[GHz] 

^Ar(vo)’^ 

<^5 

U<K] 

< 



Haslam 

Effelsberg/Jodrell/Parkes 

0.408 

1.000 

3x 10^ 




K band 

WMAP 

22.8 

0.987 

6.00 


1.007 

1.549 

30 GHz 

Planck LFI 

28.4 

0.979 

2.49 

1.013 


3.585 -0.485 

Ka band 

WMAP 

33.0 

0.972 

4.44 



... -0.788 

Q band 

WMAP 

40.6 

0.958 

3.86 



... -1.342 

44 GHz 

Planck LFI 

44.1 

0.951 

2.80 



-9.571 -2.457 

Vband 

WMAP 

60.8 

0.910 

4.32 




70 GHz 

Planck LFI 

70.4 

0.881 

2.22 




Wband 

WMAP 

93.5 

0.802 

5.14 




100 GHz 

Planck HFI 

100 

0.777 

0.86 




143 GHz 

Planck HFI 

143 

0.592 

0.37 

-1.052 

-1.088 

6.194 3.648 

217 GHz 

Planck HFI 

217 

0.334 

0.51 




353 GHz 

Planck HFI 

353 

0.075 

1.69 

0.039 

0.038 

-0.208 -0.126 

545 GHz 

Planck HFI 

545 

0.006 

0.5U 




857 GHz 

Planck HFI 

857 

6 X 10-5 

0.48" 





^ Reference frequency. 

^ Conversion factor from differential thermodynamic temperature to Rayleigh-Jeans brightness temperature. 

^ Median rms per beam for the 1° smoothed Stokes I maps. 

^ Weights for the linear combination images shown in Figs. 1(a) and 2(a) (wq), Fig. 2(b) (wi), Fig. 3(a) (W2), and Fig. 3(b) (W 3 ). 
^ Units MJysr“h 


Instrument (LFI) and High Frequency Instrument (HFI). Details 
of the map preparation, including smoothing to 1° full-width 
half-maximum (FWHM) resolution and correction of residual 
offsets and gain errors, are given in Planck Collaboration X 
(2016). The products are at HEALPix A^side = 256. In addition, 
for polarization analysis we have prepared maps at 2° resolution, 
using the same methodology as for the 1° maps. Component sep¬ 
aration depends critically on the accuracy of the input data, and 
so we briefly review data quality issues here. 

Thermal noise. This is accurately characterized in the Planck 
and WMAP data. For both missions, the noise is quite vari¬ 
able over the sky. We have calculated the thermal noise in our 
smoothed maps via propagation of errors, assuming the noise is 
independent between pixels at our full resolution. After smooth¬ 
ing, the noise in nearby pixels is strongly correlated. 

Confusion. Confusion noise denotes the uncertainties caused 
by faint compact sources within each beam. It has been esti¬ 
mated from source counts made at much higher angular reso¬ 
lution. To some extent, confusion is automatically fitted out as 
part of component separation. This would be exactly true if the 
source had a spectrum that could be well fitted by our model, and 
was non-variable. In fact, typical sources responsible for confu¬ 
sion in our maps are blazars, characterized by strong variability 
and “flat” {fi ^ -2) spectra, which are similar to free-free emis¬ 
sion, but may contain multiple peaks (e.g., Planck Collaboration 
XIII 2011). Hence a fraction of the confusion is an effective 
noise component in our fitting. 

Calibration. Uncertainties for both WMAP and Planck due 
to calibration are about 0.2% for most frequencies and this 
contributes negligibly to our error budget. The two highest- 
frequency HFI channels have raw uncertainties of < 10 % 
(Planck Collaboration VIII 2016), and corrections were derived 
as part of the component separation fit; these uncertainties have 


very little impact on the low-frequency components discussed 
in this paper. The global calibration uncertainty of the Haslam 
et al. (1982) map is about 10%, but this only impacts our es¬ 
timate for the break in the synchrotron spectrum. There is also 
an additional source of error that arises from the broad band- 
widths of the Planck and WMAP detectors. We can take these 
into account by applying “colour corrections” that depend on the 
spectral response of the instrument and the spectrum of the sky 
signal (Jarosik et al. 2011; Planck Collaboration II 2014; Planck 
Collaboration IX 2014; Planck Collaboration II 2016; Planck 
Collaboration VII 2016). For Planck LFI these corrections are 
typically 1-3 % while they are 5-10 % for HFI. The residual er¬ 
rors from such corrections are likely to be < 1 %, especially for 
the LFI channels of interest here. We implement WMAP colour 
corrections using the same approach as for the LFI, allowing for 
the secular drift of the centre frequency (Bennett et al. 2013a) 
and the small bandpass shifts derived in Planck Collaboration 
X (2016). Uncertainties due to residual beam asymmetry in the 
smoothed Planck data will be, relative to other uncer¬ 

tainties, minimal because we have used the deconvolved WMAP 
data and the additional smoothing to 1° FWHM resolution will 
symmeterize the beam. 


Offsets. Neither Planck nor WMAP measured the absolute zero 
level in total intensity, and so there is an arbitrary offset in each 
map. There is also a small but obvious residual dipole in dif¬ 
ferences between Planck and WMAP maps. In the analysis in 
Sect. 3 we use the approach of Wehus et al. (2014) that uses 
T-T plots to estimate the dipole and monopole amplitudes for 
each map (exact values were updated based on the 2015 release 
Planck maps); the approach used for the Commander analysis is 
described in Planck Collaboration X (2016). The zero levels in 
the polarization maps (Stokes Q and U) are much better defined, 
and are recovered with negligible error by Planck. In contrast, 
two large angular scale modes in the WMAP Q and U maps 
are sensitive to zero-level uncertainties in the raw data. These 
“poorly constrained modes” were down-weighted by the WMAP 
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team in their power spectrum analysis (Bennett et al. 2013a), but 
this approach is not available when attempting to model the spec¬ 
trum pixel-by-pixel. They are strong enough to seriously affect 
the absolute polarized brightness over most of the sky, except for 
the regions of strongest polarized intensity, such as the Galactic 
plane and the North Polar Spur. 

Polarization leakage. Various processes cause leakage of to¬ 
tal intensity into the polarization signal, notably mismatch be¬ 
tween the beamshapes for the two polarizations in each horn, 
and bandpass mismatch (e.g., Leahy et al. 2010). The WMAP 
scanning strategy allows a rather accurate correction for these ef¬ 
fects, but a more indirect approach is needed for Planck (Planck 
Collaboration II 2016) and the corrections applied are believed 
to be accurate to a few tenths of a percent of the total intensity. 
This is a significant problem along the Galactic plane, where the 
fractional polarization is very low and so the residual bandpass 
leakage error is relatively high. Beam mismatch leakage is also 
strong on the plane, due to the steep intensity gradients there. But 
for most of the sky, the residual leakage is essentially negligible. 

Un model led emission. This includes the Sunyaev-Zeldovich 
(SZ) effect over most of the sky,^ most spectral lines, and the 
spectrum of extragalactic point sources. Hundreds of transitions 
are known in the part of the spectrum of interest here, but over 
most of the sky they are expected to be very weak compared to 
the continuum in our broad bands. However, they are not always 
negligible: rich line spectra are generated in photo-dissociation 
regions (PDRs) surrounding H ii regions, and also in the X-ray 
dominated region around the Galactic centre (Takekawa et al. 
2014). While this spectral contamination is certainly unusual, 
there are no blind spectral line surveys in most of our bands that 
would allow us to set reliable upper limits to line contamina¬ 
tion. The Coininander fit actually includes line emission (CO and 
HCN) in four of our bands. 

2.2. Ancillary data 

Below 20 GHz, several radio maps are available. Reich et al. 
(2004) have presented a preliminary map at 1.4 GHz, but the 
combination of instrumental artefacts and calibration uncertain¬ 
ties in the survey do not meet our stringent requirements on the 
error budget, which also caused us to drop a number of Planck 
datasets as described in Planck Collaboration X (2016). We do 
use the 408-MHz map of Haslam et al. (1982), as re-processed 
by Remazeilles et al. (2015) to reduce instrumental striping and 
provide better point-source subtraction than previous versions. 
The noise in the Haslam et al. (1982) map is not characterized in 
detail, but in any case is much less than residual scanning arte¬ 
facts. We assume an rms error of 1 K, including both effects. We 
also assume a 10 % uncertainty in the overall intensity scale. 

The Galactic synchrotron spectrum shows significant curva¬ 
ture below 10 GHz, while free-free absorption sets in for bright 
HII regions below a few hundred MHz. We therefore do not con¬ 
sider maps at frequencies below 408 MHz. 

We use the full-sky dust-corrected Ha map of Dickinson 
et al. (2003) to compare with our free-free solution at high lat¬ 
itudes where dust absorption is relatively small. The fraction of 
dust lying between us and the Ha-emitting gas is assumed to be 

^ The largest SZ signals, due to the Coma and Virgo clusters, are 
explicitly fitted by Commander, but this is only feasible because other 
foreground components are very faint in these regions. 


1/3 (/d = 0.33). We also use the Finkbeiner (2003) Ha map to 
test the results. At low latitudes, we make use of the radio recom¬ 
bination line (RRL) survey recovered from the HIPASS data at 
1.4 GHz (Alves et al. 2014). This survey detects RRLs from ion¬ 
ized gas along the narrow Galactic plane with a resolution of 15'. 
Although of limited sensitivity, it provides the radio equivalent 
of Ha maps without absorption along the line of sight. 

We use the IRIS reprocessed far infrared (FIR) map from 
the IRAS survey at 100, 60, 25, and 12 p,m (Miville-Deschenes 
& Lagache 2005). The intensity scale is good to 13 % and in¬ 
strumental noise is effectively insignificant after smoothing to 
1° resolution. We also use the Hi map from the LAB survey 
(Kalberla et al. 2005). 

We also make use of the ROSAT diffuse sky survey maps. 
The original images provided by the Max-Planck-Institut fiir 
Extraterrestrische Physik (MPE)"^ are stored as six panels that 
between them cover the whole sky, with slight overlaps between 
adjacent panels. The pixel size is 12'. The 0.1-2keV energy 
range is divided into six bands. The publicly available FITS im¬ 
ages are derived from the analysis by Snowden et al. (1997) of 
the 1990-1991 ROSAT survey observations, and do not include 
the make-up data collected in 1997 and presented by Freyberg & 
Egger (1999); consequently they show narrow strips of missing 
data. After correcting erroneous field centres in the MPE EITS 
headers, we found the pixel coordinates in the original images of 
HEALPix (Gorski et al. 2005) pixels at Aside = 2048, so that the 
original pixels are highly oversampled. Nearest-neighbour inter¬ 
polation was used for placing the intensity and uncertainty val¬ 
ues into the HEALPix grid. We averaged the data and uncertain¬ 
ties onto an Aside = 256 pixel grid (pixel size 13.'7, comparable 
to the original), weighting each raw pixel by the area contributed 
to the output HEALPix pixel. We then merged the HEALPix maps 
of the six original panels into a single all-sky map, using simple 
averaging for pixels covered by two or more panels.^ 

3. A first look at the low-frequency foregrounds 
using constrained internal linear combinations 

Of all the emission mechanisms contributing to the cm- 
wavelength sky, the CMB is the most accurately measured, 
since the data are calibrated directly on the CMB dipole, and 
it is also the only component with a precisely known spec¬ 
trum. It can therefore be straightforwardly eliminated by con¬ 
structing linear combinations of the individual frequency maps. 
Therefore, as a first look, we apply a linear combination to the 
y^MAPjPlanck maps to investigate the range of spectral indices 
and morphologies at frequencies around 20-50 GHz. Our ap¬ 
proach (Appendix A) is similar to, but simpler and more limited 
than, standard ILC techniques. 

Eigure 1(a) shows the total foreground emission in the 
Planck 30 GHz band (reference frequency of 28.4 GHz), con¬ 
structed via 

^foreground = >^30^30 + >^^143^143 + (1) 

where the subscripts refer to channel frequencies, the weights 
Wi are given in Table 1 and, as always in this section, are nor¬ 
malized to return a = -3 power law unchanged. We use the 

4 www.xray.mpe.mpg.de/rosat/survey/sxrb/12/ass.html. Although 
these web pages display GIF images including the 1997 make-up 
observations, the downloadable FITS images include only the 1990-91 
data. 

^ All-sky HEALPix maps of the six ROSAT bands created as de¬ 
scribed here are available at www.Jb.man.ac.uk/research/cosmos/rosat/. 
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Fig. 1. (a) CMB-nulled map at 28.4 GHz, constructed using the ILC method (see text), where in this map, ^ = -3 emission is 
unchanged by construction; (b) 408 MHz map (Remazeilles et al. 2015), strongly dominated by synchrotron emission; (c) 545 GHz 
Planck map, strongly dominated by thermal dust emission; and (d) Her map (Dickinson et al. 2003). An asinh(/) colour scale is 
used for all images, where I is the intensity in the units indicated. The coordinate along the colour bar is linear in 1 near zero and 
becomes logarithmic (ln(2/)) at high intensity. All maps are at 1° resolution. 


Planck HFI 143 GHz band as our primary model of the CMB, 
since it has the lowest noise and is dominated by only two com¬ 
ponents, the CMB itself and thermal dust emission. The latter 
is cleaned to adequate precision for display purposes using the 
HFI 353 GHz map. The dust is modelled as a uniform mod¬ 
ified blackbody spectrum with (3^ = 1.51 and = 19.6 K 
(Planck Collaboration Int. XXII 2015). This combination frac¬ 
tionally diminishes the amplitude of free-free emission, but only 
by 0.7 %, which is negligible for display purposes. 

Figure 1 also shows our ancillary maps at 408 MHz (trac¬ 
ing synchrotron emission), Ha (tracing free-free emission), 
and the Planck 545 GHz map (tracing thermal dust emission). 
Features of all these templates appear in the 28.4 GHz fore¬ 
ground map: the brightest features along the Galactic plane are 
mostly free-free emission; at high latitudes in the northern hemi¬ 
sphere the North Polar synchrotron spur is prominent; while 
obvious dust-related features include the p Oph complex, the 
Chameleon/Musca complex of cold clouds arcing around the 
South Celestial Pole, the Polaris Flare, and the comet-shaped 
R CrA molecular cloud. In addition to the diffuse components, 
there are thousands of point sources. In fact, at this resolution 
the map is strongly confusion-limited (Fig. 2); the median ther¬ 
mal noise is 2.5 p,K, but the rms in patches a few degrees across 
is about 10 p,K, even in regions far from the plane where cat¬ 
alogued sources have been avoided. Consequently a signal is 
detected significantly in almost every pixel, and only a few SZ 


decrements (notably from the Coma Cluster) give a negative map 
intensity. 

After the CMB, the next best determined spectrum is free- 
free emission. Assuming local thermodynamic equilibrium, the 
spectral shape depends, very weakly, on a single parameter, the 
electron temperature Tq\ we use the approximation of Draine 
(2011) detailed in Sect. 4.1 below. For plausible values, 3000 < 
Te < 15000 K, the ratio of observed brightness between our 30 
and 44 GHz bands varies only between 2.52 and 2.49 (includ¬ 
ing the effect of colour corrections). We can therefore null this 
component, taking = 7500 K, by differencing CMB-corrected 
maps at 28.4 and 44.1 GHz. Figure 3(a) shows the result. 

The drawback of this procedure is that the thermal noise is 
amplified by a factor of 11, partly because the differencing in¬ 
troduces noise from the 44 GHz image, and partly due to the 
renormalization. The situation can be considerably improved by 
including WMAP data: Fig. 3(b) shows a similar analysis that 
includes the three lowest-frequency WMAP bands, where we 
maximize the S/N ratio of our reference p = -3 spectrum as well 
as eliminating CMB, free-free, and thermal dust. In this case the 
S/N ratio is 5.5 times worse than in the original 28.4 GHz map; 
the better performance is partly due to the larger frequency ratio 
(22.8/44.1 vs. 28.4/44.1), which reduces the unwanted cancel¬ 
lation of the synchrotron emission. 

Emission with a power-law spectrum steeper than free-free, 
i.e., p < -2.1, will show up as positive in the free-free nulled 
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Fig. 2. Left: (a) South polar region of the CMB-nulled 28.4 GHz map. Overplotted circles are the positions of 30-GHz sources from 
the second Planck Catalogue of Compact Sources (PCCS2, Planck Collaboration XXVI2016). The 90 % completeness limit for the 
catalogue (426 mJy) corresponds to a peak brightness of AT 50 p,K at the 1° resolution of this map. Right: (b) same region in a 
CMB-nulled version of the WMAP K-band map (22.8 GHz). The diffuse arc in the range 30° < I < 150° is part of Loop II. The 
colour scales are linear here. 


maps, and these show (albeit with higher noise) the character¬ 
istic synchrotron and AME features noted in our description of 
Fig. 1 above. Most of the extragalactic sources are nearly nulled, 
since their spectrum is similar to free-free; those remaining are 
mostly optically-thin synchrotron sources, but sources that var¬ 
ied between the WMAP and Planck mean epochs also show up 
(as negative residuals if they were brighter during the Planck 
observations). 

In this double difference image the assumed thermal dust 
spectrum begins to have a significant impact; if we use = 1.66 
and Td = 19.0 K (Planck Collaboration Int. XIV 2014) the inner 
Galactic plane becomes about 20 % brighter. Given that Td is 
known to vary significantly between the inner and outer Galaxy 
(e.g., Planck Collaboration XI 2014), more sophisticated mod¬ 
elling is required to obtain the most from the data. 

Nevertheless, the ILC approach does reveal an important 
complexity of the low-frequency foreground spectrum in a 
straightforward way, namely as prominent extended regions of 
negative intensity, notably coincident with the Sh2-27 H ii region 
at (/, b) = (6°, 24°), and in the southern part of the Gum nebula, 
around (/,/?) = (257°, -14°). These are regions with excess flux 
compared to the free-free spectrum near 40 GHz. They can be 
fitted with a spinning dust spectrum shifted to peak at higher 
frequencies than standard models, but still within the range of 
plausible interstellar parameters (e.g., Planck Collaboration XX 
2011; Planck Collaboration Int. XV 2014).The Sh2-27 spinning 
dust is likely to be associated with the famous translucent cloud 
in which interstellar molecules were first discovered, as absorp¬ 
tion lines in the spectrum of the ionizing star f Oph (Adams 
1941). In fact, many bright Hii regions show negative residuals, 
e.g., the California nebula (/, b) = (161°, -13°), Orion A and B 
(/,/?) = (210°,-19°) and (207°,-17°), and many of the bright 
nebulae along the Galactic plane. 

These residuals are not due to errors in our assumed val¬ 
ues; pushing Te as high as 15000 K makes almost no difference 
to the spectrum and so leaves all these features in place. They are 


also insensitive to changes in estimated colour corrections be¬ 
tween Planck Collaboration II (2014) and Planck Collaboration 
II (2016). They become less prominent for steeper assumed ther¬ 
mal dust spectra (Pf), but this is mainly due to the impact of 
on the underlying Galactic ridge emission noted above. 

The systematic negative residuals suggest that high- 
frequency spinning dust may often be associated with the PDRs 
around Hii regions. Dobler & Finkbeiner (2008b) also iden¬ 
tified a potential high-frequency spinning dust component in 
WMAP data that was correlated with Ha from the Warm Ionized 
Medium. An alternative explanation in at least some cases (in 
particular the Galactic centre) might be significant line emission 
in the 39-46 GHz region common to WMAP Q and the LFI 44- 
GHz bands; however well-known lines in this band, such as SiO 
(1-0) and methanol masers, do not seem to be bright enough to 
cause the effect (Jordan et al. 2015, Jordan priv. comm.). The 
implications of the wide range in AME peak frequency are dis¬ 
cussed in Sect. 4.2. 

From this preliminary analysis we can draw two further 
lessons. First, separating components with similar spectra makes 
great demands on sensitivity; even separating free-free emis¬ 
sion ip = -2.1) from synchrotron {p ^ -3) has dramatically 
reduced our effective S/N ratio. The synchrotron spectral in¬ 
dex varies by only a few tenths across the sky, and to map it 
with useful precision we need a substantial gain in sensitivity. 
This has been achieved mainly by working at low resolution 
(Fuskeland et al. 2014; Vidal et al. 2015); in the current paper 
we use a fixed template for the synchrotron spectrum (Sect. 4.3 
below). 

Second, low-frequency foregrounds have complicated spec¬ 
tra: at a minimum, we must solve for the amplitudes of free-free 
emission, synchrotron emission, and AME, and for the latter the 
peak frequency and spectral width. However, even including the 
WMAP data, we have only five frequency maps between 22 and 
44 GHz to solve for these five parameters. With no spare degrees 
of freedom, we would be susceptible to both errors in the as- 
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(a) 



-110 087 X 50 


Fig. 3. Maps with CMB, free-free, and thermal dust emis¬ 
sion nulled, scaled so that a yS = -3 power-law spectrum 
is unchanged in amplitude. Top: (a) Planck-on\y map (28.4, 
44.1, 143, 353 GHz), scaled to 28.4 GHz. Bottom: (b) combined 
Planck and WMAP. Weights used are given in Table 1 . An asinh 
colour scale is used and an offset of 80yuKcMB has been sub¬ 
tracted to enhance the contrast. 

sumed spectral model and any artefacts in the data. We therefore 
need to use maps at higher and lower frequencies. 

Above 44 GHz we have the Planck maps at 70 and 100 GHz 
and WMAP V- and W-bands (60 and 94 GHz). There is signifi¬ 
cant line contamination in the W- and 100 GHz-bands, which has 
to be included in the component separation. In all these maps the 
strongest foreground component is the low-frequency tail of the 
thermal dust emission, for which it is necessary to fit at least the 
temperature and probably the emissivity index (3^ as well. As 
a result, including these higher frequency maps gives us more 
data points but also require us to solve for more spectral pa¬ 
rameters. Moreover, since the foreground emission is weakest at 
these frequencies, the maps have relatively poor S/N ratio. The 
Coinmander approach discussed in the next section is our best 
attempt to date to handle these complexities. 

4. Total intensity foregrounds with Commander 

The analysis of Sect. 3 demonstrated the complexity of the low- 
frequency foreground spectrum, but for that very reason was un¬ 
able to separate the various diffuse foreground components that 
contribute to the maps. A more sophisticated method is required, 
and numerous algorithms exist that utilize the spectral and/or 
morphological information in the data. One of the most powerful 
algorithms is parametric fitting via Gibbs sampling, which has 



Ha [R] 

Fig. 4. Top: T-T plot of the free-free amplitude at 22.8 GHz 
against the Ha intensity [R] in the Gum Nebula (see Fig. 5 for 
the location of the region). Bottom: T-T plot for 85 % of the sky. 
Regions that show large dust absorption have been masked out, 
mostly within \b\ < 7°. 


been implemented in the Commander code (Eriksen et al. 2008). 
The details of the specific implementation to these data are de¬ 
scribed in Planck Collaboration X (2016). In this section, we 
discuss the models for the low-frequency components used by 
Commander, and compare the derived foreground products with 
previous results such as WMAP component-separation products. 
We also compare the results with expectations, either based on 
theory or extrapolations from ancillary data. 

4.1. Free-free 

The free-free radio continuum brightness can be estimated by the 
use of recombination lines. The most widely used is the Ha line 
(/1656.28 nm), which traces warm ionized gas and is proportional 
to the emission measure (EM), in the same way as free-free radio 
continuum (Dickinson et al. 2003). Thus Ha is a reliable estima¬ 
tor provided that the gas is in local thermodynamic equilibrium 
(LTE), that there is no significant dust absorption along the line- 
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of-sight, that there is no significant scattered component of Ha 
light, and that the electron temperature is known. The relation¬ 
ship between the free-free brightness [/iK] and Her intensity 
lua [R] for the optically thin limit is given by 

^ = 1512 io 0 029n y-zo ^ ( 2 ) 

where 74 is the electron temperature in units of 10"^ K, vghz the 
frequency in GHz, and gs is the Gaunt factor, which takes into 
account quantum mechanical effects. In the pre-factor we have 
included a factor of 1.08 to account for the Heii contribution 
that adds to the free-free continuum.^ Draine (2011) provides an 
approximation to gff that is accurate to within about 1 % in the 
Planck frequency range: 


gff(y, T^) = In 



5.960 - — ln(vGHz7’;^^2) 

71 ^ 



( 3 ) 


At high Galactic latitudes, the absorption by dust is a small 
effect (< 0.1 mag), and can be corrected for to first-order. The 
electron temperature is known to vary throughout the Galaxy, 
but has typical values of (7500 + 1000) K in the local diffuse 
ISM (Shaver et al. 1983; Paladini et al. 2004; Alves et al. 
2012). This corresponds to a theoretical radio-to-Ha ratio of 
(ll.l+0.9)/iKR-i at 22.8 GHz, assuming Te = (7500+1000) K. 
However, previous template fitting results have indicated that 
the measured ratio is lower than this value, at 7-10yuKR“^ 
(Davies et al. 2006; Ghosh et al. 2012). 

This apparent discrepancy can be resolved with a lower elec¬ 
tron temperature (T^ ^ 4000 K), although this is outside the 
range expected for our position in the Galaxy. Dust absorption 
cannot explain it since it has the opposite effect: it would in¬ 
crease the derived Ha emissivity resulting in an even lower ra¬ 
tio. The most popular explanation is that a fraction of the ob¬ 
served Ha intensity is from light scattered by dust grains (as¬ 
suming the gas is close to being in LTE). Previous estimates 
of this fraction suggested relatively low values (around 10 %), 
while more recent analyses indicate that as much as half of the 
high latitude Ha intensity could be scattered (Witt et al. 2010). 
Other analyses give values around the 10-20 % level (Wood & 
Reynolds 1999; Brandt & Draine 2012; Barnes et al. 2014). We 
attempt to estimate the scattered fraction in Sect. 4.1.1. We be¬ 
gin by comparing our derived free-free amplitude with the ob¬ 
served Ha. We note that there are essentially no constraints on 
Tq from the Commander solution, which is based on variations 
in the free-free spectral index, except in the inner Galactic plane 
(see Sect. 4.1.2). 

Figure 4 (top panel) shows slT-T plot of the free-free am¬ 
plitude against the Ha intensity for the Gum nebula region 
(which is known to be dominated by free-free emission), as de¬ 
fined in Fig. 5. A very good correlation can be seen, which in¬ 
dicates the robustness of the Commander free-free amplitude. 
The best-fitting slope is (8.9 ± 0.9)yuKR“^. This corresponds 
to an electron temperature of (5200 ± 900) K, which is consis¬ 
tent for the temperature measured around the same region using 
Hydrogen RRF by Woermann et al. (2000). We also measure 
the free-free-Hcr correlation using a full-sky map, masking the 
more dusty regions (Fig. 4, bottom panel). The mask was de¬ 
fined by ignoring all pixels with an optical depth at 353 GHz, 

^ We note that in Planck Collaboration X (2016), the equation for the 
free-free continuum brightness does not take into account this additional 
contribution from helium. 



Fig. 5. Commander free-free full-sky map at 22.8 GHz. The box 
shows the region we selected in the Gum nebula for the scatter 
plot shown in Fig. 4. 


T 353 , larger than 1.5 x 10“^ (this corresponds to an extinction at 
the Her wavelength of A(Hcr)=0.54 mag, assuming a reddening 
value Ry = 3.1; Planck Collaboration XI 2014), which results 
in a masked area of approximately 15 % of the sky, mostly on 
the Galactic plane. In the region studied, we find a best-fitting 
slope between the Commander free-free map and the Her map 
of (8.0 ± 0.8)yuKR“^. This is close to the value for the Gum 
nebula, and is consistent with previous analyses of the free-free- 
to-Her ratio at this frequency (Banday et al. 2003; Davies et al. 
2006; Ghosh et al. 2012). This gives us additional confidence 
that the component separation has worked relatively well for the 
free-free component. Assuming the theory is correct, this lower 
value for the free-free to Ha ratio suggests a lower value for Tq, 
of around 4500 K. However, it is likely that a fraction of the Her 
intensity is from scattering of Her light from the Galactic ridge 
by dust grains. 

4.1.1. High latitude Her scattering 

In this section we estimate the fraction, /scatt, of Her light that is 
scattered by dust grains. This is believed to be responsible for the 
low free-free to Her ratio discussed in the previous section, and 
could reduce the intrinsic Her intensity by up to 50 % (Witt 
et al. 2010). To do this, we must constrain Tg, since there is a 
one-to-one degeneracy between Tq and /scatt- For high latitude 
sight-lines away from the Galactic plane, we can assume that 
the bulk of the gas is at the local value Tq = (7500 ± 1000) K, 
thus breaking the degeneracy. We initially ignore variations in 
dust absorption, since we have masked out the regions that are 
most affected. We test the effect of this assumption by making a 
nominal correction for dust absorption. 

First, we use the theoretical proportionality constant that re¬ 
lates the Her intensity with the Commander free-free emission, 
(11.1 ± 0.9)/iKR“^ to scale the free-free template to Rayleigh 
units. We note that the uncertainty accounts for a +1000K un¬ 
certainty in Tq. We subtract this template from the Her map to 
obtain the Her emission that does not come from in-situ recom¬ 
bination, i.e., originating from Her light scattered by dust grains. 
The top panel in Fig. 6 shows the residual map, with the pixels 
masked out shown in grey and the colour scale being truncated to 
±2 R, which saturates the brightest regions closer to the Galactic 
plane and the Gould Belt. A feature of this residual map is the 
blue central region, which resembles the microwave haze noted 
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- 2.0 2.0 R 



Fig. 6. Top: residual Ha emission after subtraction of the scaled 
free-free component. Bottom: scatter plot of the residual Ha in¬ 
tensity against the thermal dust optical depth at 353 GHz out¬ 
side the mask shown above, corresponding to A(Ha) < 0.53 
(^353 < 1 X 10“^). The best-fitting linear relation is shown as 
a dashed line. 


by Finkbeiner (2004), Dobler & Finkbeiner (2008a) and Planck 
Collaboration Int. IX (2013). This is not unexpected, since the 
haze was discovered by assuming that Ha is an accurate free- 
free template. This is an indication that component separation 
is ambiguous in this complex region, since mixed emission from 
the Gould Belt system makes the separation between AME, free- 
free, and synchrotron particularly difficult. There may also be 
component separation issues in other regions of the sky at a 
lower level, which could have a significant impact on the these 
residuals. 

However, we can still set an upper limit for the amount of 
scattered Ha light, based on the value of the free-free to Ha 
ratio measured over the sky. The mean value at 22.8 GHz that 
we measure over most of the sky is (8.0 + 0.8)yuKR“^ From 
this, if we assume that the difference relative to the theoretical 
value of (11.1 + 0.9)yuKR“^ is due to an excess of Ha emis¬ 
sion from scattering, we can estimate the scattering fraction to 
be /scatt = (28 + 12) %. This is an average value for the high lat¬ 
itude sky, so individual regions with difiTerent electron tempera¬ 
ture might be present a level higher than this. It is consistent with 
previous estimates of around 20 % of scattered light on average 


at high latitudes (Wood & Reynolds 1999; Witt et al. 2010; Dong 
& Draine 2011; Brandt & Draine 2012; Barnes et al. 2014). If we 
repeat the calculation using a dust-corrected Ha map, assum¬ 
ing that 1/3 of the dust lies in front of the Ha-emitting gas (see 
Dickinson et al. 2003), we find /scatt = (36 + 12) %. 

Our constraints are clearly not very tight, and are dependent 
on our assumption on the mean electron temperature and on the 
level of dust absorption. Given the relatively large uncertainties, 
an electron temperature of about 5000 K would be enough to 
bring the ratios within the Icr uncertainty of the observed range. 

We now explore the correlation between the residual Ha 
emission and the dust optical depth at 353 GHz. Scattering by 
dust grains should produce a positive correlation between these 
two maps. We plot this in the bottom panel of Fig. 6 , where for 
low values of the opacity, T 353 < 1 x 10 “^, the correlation is in¬ 
deed positive. The large dispersion of the points in the figure is a 
consequence of the noise and variations of electron temperature 
at high Galactic latitude. In spite of that, a significant correlation 
is observed. This suggests that the effect of scattered Ha light 
might be real and that its effects must be taken into account to 
determine precise values of T^. But we remind the reader that 
this could also partly be explained by the regions corresponding 
to 0.5-1 R (at high Galactic latitudes, \b\ > 40°, where the pos¬ 
itive correlation is observed) having a lower than average value 
for Te. 

We do not find strong evidence for high latitude regions with 
a scattered light fraction of around 50 %. We also do not see an 
anti-correlation of the free-free-to-Ha ratio with the thermal dust 
optical depth (a proxy for column density of dust), as would be 
expected if there were a significant fraction of scattered light. We 
hypothesize that the electron temperature at our exact Galactic 
position could be slightly lower than expected (say 6000 K), 
which would reduce /scatt to ~ 10 %. 

4.1.2. Free-free emission in the plane 

Along the Galactic plane, the Ha line is highly absorbed by dust 
and thus fails to provide a reliable measure of the free-free emis¬ 
sion (Dickinson et al. 2003). In this region of the sky, we can 
use RRLs to estimate the thermal continuum emission and com¬ 
pare with the Commander free-free solution. The integrated RRL 
brightness is proportional to the free-free continuum in exactly 
the same way as for the optical Ha recombination line. 

We use the data from the Hi Parkes All-Sky Survey 
(HIPASS; Barnes et al. 2001) of the southern sky, which have 
been re-analysed to extract RRLs at 1.4 GHz to produce the 
first continguous RRL survey of the Galactic plane (Alves et al. 
2014). The full survey covers the longitude range 196° ^ 0° ^ 
52° and \h\ < 5°, with an angular resolution of 14.'4. The inten¬ 
sity calibration is good to better than 10 %. 

Figure 7 shows maps of the free-free emission from the 
Commander fit and from the RRLs. The free-free amplitude 
is converted from the RRL temperature assuming a simple 
model for Tg, based on an assumed temperature gradient of 
(500 + 100)Kkpc“^ with distance from the Galactic centre 
(Alves et al. 2014). The Galactocentric distance at each longi¬ 
tude is estimated based on the Fich et al. (1989) rotation curve 
and using the central velocity of the RRLs. Obviously this is 
a crude approximation, since the RRL is a line-of-sight inte¬ 
gral weighted by and individual lines of sight could have a 
range of values. We estimate that the electron temperatures for 
a given line-of-sight are uncertain at the + 1000 K level, which 
corresponds to an uncertainty in the free-free amplitude of up 
to 20 %. Thus the overall uncertainty on the free-free amplitude 
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Fig. 7. Maps of the free-free emission along the Galactic plane, 
separated into four different longitude ranges. For each /-range, 
the RRL and Commander estimates are shown in the top and bot¬ 
tom panels, respectively. The maps are in units of K at 1.4 GHz 
and at 1° resolution. There is a single contour for each pair of 
panels, set at a value corresponding to the minimum tempera¬ 
ture of the 2 % brightest pixels of the RRL map, corresponding 
to 7.4, 9.6, 4.4, and 0.4 K, from top to bottom, respectively. The 
white pixels at the location of some bright Hii regions in the 
maps of the first, third and fifth rows, correspond to saturated 
pixels in the RRL survey. 


from the RRL data is 25 % at most. Nevertheless, there is an 
overall good level of agreement between the two estimates, both 
in terms of morphology and amplitude. A combination of bright 
compact sources (H ii regions) and diffuse emission within about 
+ 1° of the plane contribute, particularly in the inner Galaxy 


(/ = 300° ^ 0° ^ 60°). In the outer Galaxy (/ < 300°) the 
free-free emission is much weaker and only a few moderately 
bright H ii regions can be seen. 

The scatter plots of Fig. 8 indicate that the two data sets 
are highly correlated, which suggests that the component sep¬ 
aration of the free-free emission from the other components has 
been relatively successful. However, there are clearly changes in 
slope from region to region. In general, the slopes are greater 
than I, indicating that the Commander free-free map is, on aver¬ 
age, brighter that the predicted RRL map. The best-fitting slopes 
are typically 50 % higher than the predicted values. 

Figure 9 shows longitudinal slices across the Galactic plane 
(b = 0°) at 1° resolution. The bright Hit regions are clearly 
visible as bright peaks. Virtually all the bright peaks in the 
Commander free-free solution have a counterpart in the RRL 
prediction. These peaks sit on diffuse emission that should be 
subtracted before a comparison of peak temperatures is made; 
large-scale (> 4°) emission will not be reliable in the RRL data 
(Alves et al. 2014). Many of the peaks have comparable ampli¬ 
tudes, while some of them are discrepant by up to a factor of 
2. We note that the Galactic centre (/ = 0°) is low in the RRL 
map due to receiver saturation artefacts. The ratio of free-free 
to RRL brightness for the brightest 6 peaks in the first region 
(/ = 358° ^ 0° ^ 52°, excluding the I = 0° peak) is 1.14 + 0.04, 
while the next 10 brightest peaks have a ratio of 1.36 ± 0.08 (the 
uncertainty is the rms scatter of the ratios). 

The discrepancy between the two free-free estimates on 
bright Hii regions, which has been noted before (Alves et al. 
2010; Planck Collaboration Int. XIV 2014; Planck Collaboration 
Int. XXIII 2015), is presumably due to a combination of ef¬ 
fects: (i) differences in the electron temperature used to estimate 
the free-free brightness; (ii) beam effects in the RRL map. The 
RRL data are calibrated in the full beam scale, such that the flux 
density of compact sources is underestimated by 10-30% de¬ 
pending on their size, although the longitude profiles of Fig. 9 
suggest that the discrepancy between the Commander and RRL 
maps across individual sources does not depend on their ex¬ 
tent; and (hi) residual emission in the Commander map, e.g., at 
(/, b) = (46?8, -0?3) associated with the HC30 supernova rem¬ 
nant. 

Compared to the brightest sources in each region, the 
weaker diffuse emission appears to have an even steeper slope, 
approaching 2, i.e., the discrepancy between the RRL and 
Commander solution is worse for very extended emission. This 
is likely to be due to excess synchrotron radiation at higher fre¬ 
quencies that has been accounted for in the simple Commander 
model by the free-free component. In particular, the synchrotron 
component is effectively modelled as a power law in frequency 
with a fixed spectral index, while there is considerable evi¬ 
dence for spectral flattening of synchrotron emission at low lat¬ 
itudes (Kogut et al. 2007; de Oliveira-Costa et al. 2008; Gold 
et al. 2011; Planck Collaboration Int. XXIII 2015). This will 
contribute to some level of excess emission at frequencies 20- 
100 GHz, which can result in an apparent increase in the free- 
free and/or AME amplitude in the Commander fits, an effect that 
can be seen as a broad background (zero-level) in the longitude 
plots of Fig. 9. We also note that the RRL data will not reliably 
trace large-scale (> 4°) emission, since the observations were 
made in 8° x 8° patches of sky. 

In summary, the Commander free-free map appears to be a 
reliable tracer of the brightest Hii regions, with an accuracy 
of around 20%. Weaker emission regions appears to be over¬ 
estimated by up to a factor of 2, relative to estimates using RRL 
data. As mentioned in Planck Collaboration X (2016), including 
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Fig. 8 . Scatter plots of the Commander versus the RRL free-free map for three regions of the Galactic plane (the first three regions 
shown in Fig. 7. The slope values given in each panel result from a linear fit to all the data points. 



Longitude [degree] Longitude [degree] Longitude [degree] 

Fig. 9. Longitude profiles of the Commander (black) and RRL (red) free-free maps, averaged within \b\ < 1°, for three regions of the 
Galactic plane (the first three regions shown in Fig. 7). The comparison is for the bright compact regions that show as peaks on top 
of a broad baseline level due to extended emission. 


new data in the frequency range 2-20 GHz will yield improved 
component separation products, particularly for synchrotron and 
AME, which in turn will improve the accuracy of the free-free 
solution. 

4.2. Anomalous microwave emission 

Commander models the AME as spinning dust, which for uni¬ 
form conditions (grain size distribution, ambient radiation) gives 
a relatively sharply peaked spectrum; specifically it uses the 
SpDust2 code (Ali-Hai‘moud et al. 2009; Silsbee et al. 2011) to 
calculate a template spectrum using parameters typical for the 
diffuse cold neutral medium (Draine & Lazarian 1998). To toler¬ 
able accuracy for present purposes, models for other phases can 
be approximated by scaling the template in frequency. 

As shown in Sect. 3, the AME spectrum is quite vari¬ 
able, with an apparent tendency to peak at substantially higher- 
than-average frequencies around some H ii regions. Low-latitude 
lines-of-sight will therefore contain AME with a range of peak 
frequencies, which will give a broader spectrum that cannot be 
fitted by a frequency-shifted template. In fact, a superposition 
of AME spectra with a distribution of peak frequencies could 
approximate a power law over the frequency range of primary 
interest here (20-70 GHz). Such a general model could not be 
constrained by the available data, but we take a step in this di¬ 
rection by using two AME components in the Commander fit, 
each based on our SpDust2 template. One component (which 
generally dominates) is scaled in frequency at each pixel, with 
a Gaussian prior for the frequency of the intensity peak, Vpi, 
of (19+ 3) GHz. The peak frequency of the other component 
is assumed to be constant over the whole sky with a value to 
be determined globally; in this Commander solution the result 


is yp 2 = 33.35 GHz. This second model component is there to 
account for the broadened spectrum of the overall AME emis¬ 
sion, and should not be interpreted as a physically distinct AME 
component; in this paper (as in Planck Collaboration X 2016) we 
only consider the combination of the two components, which we 
evaluate at 22.8 GHz to allow comparison with earlier studies. 

The all-sky map of AME as fitted by Commander is shown in 
Fig. 10. The foreground emission detected at frequencies around 
20-60 GHz in the WMAP/P/a^c^ maps appears to be domi¬ 
nated by AME (e.g., Davies et al. 2006; Ghosh et al. 2012, 
Sect. 3 of this paper). It has been detected with high signif¬ 
icance over a large area of the sky. The strongest emission is 
observed in the Galactic plane, particularly in the inner Galaxy 
(/ = 300° ^ 0° ^ 60°). The AME in the inner plane was iden¬ 
tified and discussed by Planck Collaboration XXI (2011) and 
Planck Collaboration Int. XXIII (2015). The map also shows 
clear, strong AME in many of the molecular cloud regions that 
are well known to exhibit AME, including the Perseus and p 
Ophiuchus clouds. While previous Planck AME papers (Planck 
Collaboration XX 2011 ; Planck Collaboration Int. XV 2014) fo¬ 
cused on compact AME regions, here we also look at diffuse 
AME regions away from the Galactic plane (\h\ > 10°). We 
identify strong diffuse AME in the Chameleon/Musca regions as 
well as the Orion complex, which is discussed later in this sec¬ 
tion. These regions are less likely to be affected by foreground 
degeneracies than others: regions such as the Auriga dust feature 
above the Galactic centre and a large structure around the North 
Celestial pole can also be seen in the AME map, however they 
overlap with significant sources of synchrotron emission. 

Overall there is a high level of correlation of AME with the 
thermal dust maps made at higher frequencies of Planck and in 
the infrared; however, the AME clearly cannot be accounted for 
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Fig. 10. All-sky map of AME from Commander at 22.8 GHz plotted using a Mollweide projection, with a 30° graticule and an asinh 
colour scheme. Seven regions of diffuse AME have been highlighted, and are discussed in the text. 



Fig. 11. Mask used for the AME T-T plots. The regions in grey 
are masked in the T-T plots. The region in light blue shows 
where the AME at 22.8 GHz is brighter than 1 mK; this is also 
masked when determining the best-fitting emissivity. We also 
mask the ecliptic plane (\I3\ < 20°) when analysing the IRAS and 
WISE data. 


by the Rayleigh-Jeans tail of the thermal dust emission, which 
is at least a factor of 30 below the AME component at 30 GHz. 
At low Galactic latitudes (\b\ <10°) and in bright compact and 
diffuse AME regions, the complexity of separating AME from 
CMB, synchrotron, and free-free emission is the dominant un¬ 
certainty for the Commander solution. At high latitude, where the 


AME is weak, the Commander solution is dominated by instru¬ 
mental noise. The Commander AME map also contains residual 
extragalactic source contributions when their spectrum is inter¬ 
mediate between the assumed synchrotron and free-free models 
(e.g., Centaurus A is clearly visible, as are bright radio sources 
such as 3C 273); to remove these we mask sources in the PCCS2 
catalogue at 28.4 GHz that are brighter than 1 Jy. 

We evaluate the peak frequency of the combined AME 
spectrum for each pixel. The mean and median across the 
whole sky are 20.6 and 20.4 GHz, respectively, and at \b\ > 10° 
they are 20.5 and 20.2 GHz, respectively. Apart from bright 
regions, the peak frequency map is very noisy, largely con¬ 
strained by our prior on Vpi around 19 GHz (see section 4.1 
and figure 16 of Planck Collaboration X 2016) and contains 
residual signal from point sources. Regions such as Perseus, 
p Oph and Orion have higher than average peak frequen¬ 
cies, around 25-30 GHz, in agreement with previous analy¬ 
ses (e.g., Planck Collaboration XX 2011). In some places the 
fits are dominated by the second AME component (no¬ 
tably the Sh2-27 and southern Gum nebulae). This is con¬ 
sistent with the ILC analysis of Sect. 3, although in re¬ 
ality some of these regions have even higher peak fre¬ 
quency, e.g., Planck Collaboration Int. XV (2014) found a best¬ 
fitting value of Vp of 50 GHz for the California Nebula. 
Sh2-27 has an exceptionally high the Commander fit 

(Planck Collaboration X 2016, figure 22), which we found can 
be alleviated by allowing a similar peak frequency (recall that 
in our baseline Commander model yp 2 is fitted globally to keep 
the number of free parameters acceptable). Dobler & Einkbeiner 
(2008b) previously found a bump at ~50GHz in the SED of 
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Fig. 12. T-T plots at A^side = 64, using the mask shown in Fig. 11. The best-fitting lines and emissivities shown are determined 
where the AME amplitude is less than 1 mK. The Spearman’s rank correlation coefficient (s-rank) and the Pearson’s correlation 
coefficient (r) are also shown. Top row: AME (p,K at 22.8 GHz) compared to: (a) Planck 545 GHz (MJy sr“^); (b) the Commander 
thermal dust solution at 353 GHz (jliK); and (c) T 353 (rescaled by 10^). Middle row: (d) 100 jam (MJy sr“^); (e) Hi (Kkms“^); and 
(f) CO J=2 1 (Kkm s-i). Bottom row: (g) IRAS 12 ^tm (MJy sr-^) divided by Go; (h) IRAS 25 ^tm (MJy sr-^) divided by Go; and 

(i) dust radiance calculated from the Commander products. 


Hcr-correlated emission from the Warm Ionized Medium that 
they speculated could be a high-frequency spinning dust com¬ 
ponent, which is likely the same emission we are seeing here. 
We find that in bright free-free regions (primarily on the plane, 
but also bright high-latitude regions such as M42 and the 
California nebula) the AME peak frequency correlates with the 
free-free amplitude (at 22.8 GHz) with a best-fitting slope of 
(0.066 + 0.012) GHzmK“\ which may be an indication of the 
importance of the interstellar radiation field in the excitation of 
AME carriers (e.g., Tibbs et al. 2013). 

4.2.1. All-sky correlations 

In order to look at the correlation of the AME with other com¬ 
ponents, we have degraded the AME map and other tracer maps 
to A^side = 64 (pixel widths of 55'). We define a mask, shown in 
Fig. 11, to exclude pixels where the S/N for AME is less than 
3.0 and point sources brighter than 1 Jy, except at \b\ < 5° to 


avoid masking Galactic plane sources. However, we mask five 
Galactic plane regions containing the brightest Hii complexes 
(e.g., Cygnus-X) in the sky. We use a 3° radius disc centred at 
(/,Z?)=(49?6,-0?6), (79?2,0?4), (268?2,-l?3), (287?5,-l?l) and 
(291?8,-0?9). These are particularly bright in the IRAS 12 and 
25 p,m maps, where there could be contamination by starlight 
and/or line emission. We show the T-T plots for the entire sky 
outside of the mask for different potential tracers in Fig. 12. 
We also test the effects of additionally masking the Galactic 
plane {b < 10°). We calculate the emissivity, i.e., the ratio of 
AME emission compared to AME tracers, from the slope of 
the best-fitting line as shown in the T-T plots for pixels where 
the AME amplitude at 22.8 GHz is less than 1 mK (shown in 
light blue in Fig. 11). The difference between the Commander 
maximum-likelihood and mean solutions for the AME com¬ 
ponent is around 5-6%, depending on frequency, due to the 
non-Gaussian shape of the probability distribution. We there¬ 
fore assume a conservative 10% modelling uncertainty on top 
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of the Commander component uncertainty maps; we also as¬ 
sume a 13.5 % calibration uncertainty for the IRAS 100 p,m map 
(Miville-Deschenes & Lagache 2005). 

We find an excellent correlation (r = 0.98) with the 
Planck map at 545 GHz, shown in Fig. 12a, which is pre¬ 
dominantly thermal dust emission. A mean emissivity of 
(65 + 7)HK(MJysr')-' is found ((70 + 7) |aK(MJy 
when the Galactic plane is masked), which is consistent 
over a wide range of amplitudes. The overall emissiv¬ 
ity against IRAS 100 p,m that we find when looking at 
the entire sky is (22 + 2) p,K(MJy sr“^)“^ (Fig. 12d), with 
r = 0.91; this is somewhat lower than the weighted av¬ 
erage found by Planck Collaboration Int. XV (2014) of 
(32 + 4) p,K(MJy sr“^)“\ although the latter is for a specific 
set of 27 low-latitude clouds selected for significant emission, 
rather than an all-sky average. Davies et al. (2006) used template 
fitting to estimate the emissivity of AME compared with thermal 
dust; they calculated the emissivity in terms of the “FDS model 
8 ” dust map at 94 GHz (Finkbeiner et al. 1999), which can be 
multiplied by a factor of 3.3 to obtain the equivalent 100 p,m 
emissivity (on average). Our analysis agrees well their results, 
where they found (21.8 + 1.0) p,K(MJy sr“^)“^ across the 
whole sky outside of the WMAP Kp2 mask, and an average of 
(25.7 + 1.3) p,K(MJy sr“^)"^ in the 15 regions they examined. 
We note that this good agreement in the overall high-latitude 
emissivity comes from different component-separation methods. 

For thermal emission tracers there is some curvature in the 
relation at the brightest pixels, where there is less AME for a 
given amount of thermal dust emission. This is more notice¬ 
able in the 100 p,m correlation than the 545 GHz correlation, be¬ 
cause the 100 p,m amplitude is much more affected by the tem¬ 
perature of the thermal dust grains than the 545 GHz emission 
(Tibbs et al. 2012), which increases towards the Galactic plane 
(e.g., Planck Collaboration XIX 2011). This effect is not seen in 
the T 353 correlation, however there is a much larger scatter in 
the T353 correlation, as well as a trend towards there being more 
AME per unit T353 at high optical depth. 

There is a good correlation (r = 0.84) with H i emission inte¬ 
grated in the velocity range [-450, -h400] kms“^ from the LAB 
survey (Kalberla et al. 2005), shown in Fig. 12e, with a slope of 
(210 + 20) ((280 + 30) when the Galactic plane 

is masked); there is also a noticeable upturn and increase in the 
scatter of points in this relation above an Hi integrated veloc¬ 
ity intensity of around 3 K, in agreement with previous results 
(e.g., Lagache 2003). Planck Collaboration Int. XVII (2014) 
compared Planck and WMAP data to Hi data, and found 
an emissivity of (1.75 + 0.14) x 10^ they also con¬ 

sidered AME models, which gave slightly lower emissivities 
with T353 of 1.3 X 10^ p,K. We find an emissivity against 
T353 of (8.7 + 0.8) X 10^ ^iK (Fig. 12c, (1.1 + 0.1) x 10^ ^iK 
excluding the Galactic plane), which is lower than the val¬ 
ues found by Planck Collaboration Int. XVII (2014) by a fac¬ 
tor of 1.5-2.1. However, the emissivities have not been cal¬ 
culated on the same areas of sky, and this difference is 
within the variations of E{B - y)/A^H reported in Planck 
Collaboration XI (2014). We also compare AME at 22.8 GHz 
with the Coinmander thermal dust solution at 353 GHz and 
find a slope from the regions of (0.9+0.1)KK“^ (Fig. 12b, 


^ Planck Collaboration Int. XVII (2014) give a value of 
(17 + 1) pK(10^^Hcm"^)"^ at 23 GHz, which can be converted 
to qKT 353 using E{B - y)/VH = (1.44 + 0.02) X 10“^^magcm^ 
and E{B - y)/T 353 = (1.49 + 0.03) x 10"^ mag"^ cm"^ 

(Planck Collaboration XI 2014). 


r = 0.98, (L0+0.1)KK ^ excluding the Galactic plane); this is 
also slightly lower than the value of (1.12 + 0.03) KK“^ found 
by Planck Collaboration Int. XXII (2015). 

The upturn in the Hi correlation implies that the AME is 
related to another phase that is not traced by Hi emission. 
We also look at CO emission, which will trace denser re¬ 
gions, and is a marker for Hii emission. The correlation with 
the Commander CO J=2 1 map is rather more complicated 

than the other correlations, however. The best-fitting slope is 
(180 + 20) p,K(Kkms“^)“^ (Fig. 12f), with r = 0.91, reducing to 
(150 + 20) p,K(Kkms“^)“^ when the Galactic plane is masked. 
However, at low amplitudes the emissivity is clearly a lot lower 
than the best-fitting value, and at high amplitudes it is clearly 
higher than the best-fitting value. This could be due to compo¬ 
nent separation issues; alternatively, since the curvature largely 
takes place at AME amplitudes above 1 mK, it may indicate that 
there is less AME in the densest parts of the Galaxy than at 
higher latitudes. 

We have also looked at the correlations with the IRAS 12 p,m 
and 25 p,m data, which trace the very small dust grain popula¬ 
tion that could produce AME, as well as containing PAH line 
emission. Indications of an improved correlation with these data 
have been seen on smaller scales by e.g., Casassus et al. (2006). 
However, these data are considerably contaminated by zodiacal 
light and other systematics; as such, we have masked the eclip¬ 
tic plane (|yS| < 20°) before fitting these data. The resulting AME 
emissivity against 12 p,m emission is (460+60) p,K (MJy sr“^)“\ 
with Pearson correlation coefficient r = 0.93 and Spearman’s 
rank correlation coeffient 5 '= 0.71, and against 25 p,m emis¬ 
sion is (370 ± 40) p,K(MJy sr“^)“^ with r = 0.87 and 5 -= 0.52. 
However, these will still include other systematic effects, such 
as stellar contamination and high-latitude zodiacal light. We 
nomalize the maps by the interstellar radiation field intensity, 
which we estimate by Go = (T^/17.5)^'^^ (Ysard et al. 2010), 
so that we are correlating against column density. Dividing 
by Go improves the correlation as measured by Spearman’s 
rank, which includes the curvature in the correlation, while the 
Pearson’s rank, which relies on a linear correlation, reduces. 
This improvement in correlation is as expected by the spin¬ 
ning dust model (Ysard et al. 2010 ), however the value of Go, 
and hence the correlation coefficient, is highly dependent on the 
dust temperatures and spectral indices. The resulting emissivity 
against 12 p,m is (290 ± 30) p,K (MJy sr“^)“^ (Pig. 12g, r = 0.78, 
5 '= 0.87) and against 25 p,m it is (200 ± 20) p,K(MJy sr“^)“^ 
(Pig. 12h, r = 0.77, 5 - = 0.83). We return to this topic in the dis¬ 
cussion of the Musca region, below. 

In a recent analysis, Hensley & Draine (2015) have com¬ 
pared the Commander AME maps with different potential tracers 
of AME, and find a good correlation against dust radiance. This 
is the integrated intensity from thermal dust, and it characterises 
the energy absorbed and emitted by the thermal dust particles. 
We have also compared the AME map with dust radiance (cal¬ 
culated using the Commander products and equation 10 in Planck 
Collaboration XI 2014) in Pig. 12i. We find a strong correlation 
at low amplitudes, with an emissivity against AME at 22.8 GHz 
of (5.6 ± 0.6) X 10^ p,K, and r = 0.90. At higher amplitudes (i.e., 
in the Galactic plane), the dust radiance calculated using the 
Commander products turns over, in the same way as the corre¬ 
lation with the IRAS 100 p,m data does, implying that this is due 
to temperature effects. If we use the radiance map from Planck 
Collaboration XI (2014), which includes the 100 p,m map, in the 
analysis, this effect reduces in magnitude, but is still present. 

Hensley & Draine (2015) also compare the variations in 
AME emissivities (defined by dividing the AME map by the 
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average emissivity against dust radiance) with variations in the 
fraction of the WISE 12 jLim data (Meisner & Finkbeiner 2014) 
that can be attributed to polycyclic aromatic hydrocarbon (PAH) 
emission, in contradiction with expectations from the spinning 
dust model where PAHs are though to be the spinning molecules. 
Using the Commander thermal dust amplitude and and (3^ 
maps from (Planck Collaboration XI 2014), they find no cor¬ 
relation between AME emissivities and PAH emission. We re¬ 
peat their analysis, using the WISE 12 p,m data and masking the 
Galactic plane {\b\ < 5°) and the Ecliptic plane (\j3\ < 20°) in ad¬ 
dition to the mask shown in Fig. 11. We find similar results when 
using the dust radiance maps calculated using the products from 
(Planck Collaboration XI 2014), with correlation coefficients of 
0.10 (Planck Collaboration XI2014 radiance map) and 0.23 (us¬ 
ing the Commander dust amplitude). However, if we calculate 
dust radiance using only the Commander products, then we find 
a correlation with a slope of 1.1 + 0 . 1 , with a correlation coeffi¬ 
cient of 0.52.^ Although the Commander dust radiance map will 
be biased low due to the absence of data points tracing warmer 
dust temperatures, this demonstrates the dependence of this re¬ 
sult on the quality of the radiance map; it will also depend on 
the quality of the AME map. One possibility is that this correla¬ 
tion is affected by the large angular scale structure in the Galaxy; 
as such we return to this in the discussion of the Musca region, 
below. 

In conclusion, we find the best correlation with AME at all 
amplitudes is from the 545 GHz Planck map, followed by the 
Commander dust solution at 353 GHz and the optical depth, T 353 . 
The dust radiance has a tight correlation with AME away the 
Galactic plane, but has a worse correlation in the plane. The cor¬ 
relation with 100 |Ltm is significantly affected by temperature ef¬ 
fects, which the choice of lower frequency thermal dust maps 
avoids. We find a reasonable correlation with Hi emission, al¬ 
though this does not appear to correlate with the brightest AME 
emission, and with CO emission, although this is not well-fitted 
with a single emissivity. The correlation with dust radiance is 
very good at low amplitudes, but using a single emissivity would 
over-predict the amount of AME present in the Galactic plane. 
We caution that these correlations will depend on the choice of 
mask due to large-scale biases, contamination from other emis¬ 
sion mechanisms (e.g., point sources, zodiacal light, etc.), and 
variations across the sky. 

4.2.2. Diffuse AME regions 

We now move on to consider individual diffuse regions of AME. 
We have labelled seven regions in the all-sky AME map of 
Fig. 10 that demonstrate diffuse AME away from the Galactic 
plane and are in areas with relatively high S/N ratios and clean 
component separation outputs. These include the well-known re¬ 
gions of Perseus and p Ophiuchus, as well as five new regions. In 
Fig. 13 we show the Commander solutions for synchrotron, free- 
free, AME and thermal dust, along with the CMB-subtracted 
28.4 GHz Planck data and Ha data, for the new regions. In each 
of these maps, an area has been defined to analyse the properties 
of the region using T-T plots, and sources brighter than 5 Jy in 


^ The first preprint version of Hensley & Draine (2015) used a full- 
sky WISE map that later turned out to be contaminated by Planck 
857 GHz data on scales larger than 2°, which might also affect our anal¬ 
ysis of their results. A revised version, Hensley et al. (2015), used an 
alternative approach to process the WISE data to avoid this issue; they 
found that their conclusions were unchanged. As such, we have not re¬ 
vised our conclusion here. 


Table 2. The top section presents emissivities of AME at 
22.8 GHz (|LtK) relative to the Commander thermal dust ampli¬ 
tude at 545 GHz (MJy sr“^), the IRAS map at 100 jam (MJy sr“^), 
and the optical depth at 353 GHz, T353, for the regions shown in 
Figs. 10 and 13. We also include emissivities for the EMC and 
SMC, discussed in Sect. 4.5. The bottom part presents emissivi¬ 
ties from Davies et al. (2006, D06; whole sky and region mean) 
and Planck Collaboration Int. XV (2014, XV; Perseus, p Oph; 
and the unweighted region mean) for comparison. 


Region. 

AME /545 GHz 

AME /100 ^m 

AME/T 353 


[|^K(MJysr-‘)-‘] 

[|^K(MJysr-‘)-‘; 

1 [vtK10-«] 

Rl: Perseus. 

24+7 

12.3 + 

1.9 

1.5+ 0.9 

R2: Plume . 

47+6 

18 + 

2 

7.7 + 1.0 

R3:RCrA . 

36 + 14 

50 + 12 

4.1 + 1.8 

R4: p Oph. 

40+9 

4.6 + 

0.9 

2.2+ 1.2 

R5: Musca . 

59+8 

26 + 

3 

6.9 + 1.0 

Chamaeleon . 

74+8 

22 + 

2 

11 +1.1 

R6: Orion. 

47+5 

20 + 

2 

4.7 + 0.6 

R7:/lOrionis ... 

104+11 

25 + 

3 

15 + 1.8 

EMC . 

56+6 

8.5 + 

1.0 

7.9 + 0.9 

SMC . 

30+9 

4.7 + 

1.5 

3.3 + 1.0 

Entire sky. 

65+7 

22 + 

2 

8.3+ 0.8 

|Z?| > 10°. 

70+7 

21 + 

2 

9.7 + 1.0 

XV: Perseus .... 


24 + 

4 


XV: p Oph . 


8.3 + 

1.1 


XV: Mean. 


32 + 

4 


D06: Kp2 mask . . 


21.8 + 

1.0 


D06: Region mean 


25.7 + 

1.3 



the PCCS2 catalogue have been masked out to a radius of 60', 
with some exceptions noted below. We show T-T plots between 
AME and thermal dust for three regions in Fig. 14. We use the 
best-fitting slope from these T-T plots to determine the emissiv¬ 
ity of the AME component compared to the other components. 
Results for all regions are given in Table 2. 

1. Perseus. This molecular cloud is a well-known source of 
AME (Watson et al. 2005; Planck Collaboration XX 2011). For 
the correlation analysis, we focus on an 8 ° x 8 ° patch cen¬ 
tred on (/,/?) = (160?26,-18?62). Perseus is included in 
the PCCS2 catalogue, so we unmask the point sources when 
calculating the emissivities in this region. We find a lower 
emissivity for Perseus in the Commander products than in 
Planck Collaboration Int. XV (2014); this is due to leakage of 
around 40 % of the AME emission to the free-free component, 
which also causes a decrement in the synchrotron map (see 
Fig. 33 below). 

2. Pegasus plume. An example where the AME amplitude 
traces a filament-like structure at the edge of the Pegasus con¬ 
stellation at (/, b) = (92°, -37°) and approximately 10° in length 
is shown in the first column of Fig. 13 (Region 2 in Fig. 10). 
We look at a 15° x 15° patch centred on this position. The 
plume was included in the CO catalogue of Magnani et al. 
(1985) as MBM 53-55 (G92?97-32?15, G92?97-37?54, and 
G89?19-40?94, respectively), and MBM 55 is also coincident 
with the HII region S 122 (G89? 18-41? 13) in Sharpless (1959). 
Fukui et al. (2014) have compared the CO emission to thermal 
dust emission as seen by Planck at 353-857 GHz. The struc¬ 
ture can be seen clearly in the CMB-subtracted Planck 28.4 GHz 
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Fig. 13. Commander component maps for diffuse AME regions. From top to bottom: synchrotron amplitude (mK at 22.8 GHz), 
Her (Dickinson et al. 2003, Rayleighs), free-free amplitude (mK at 22.8 GHz), AME amplitude (mK at 22.8 GHz), thermal dust 
amplitude (MJy sr“^ at 545 GHz), CMB-subtracted Planck 28.4 GHz map (mK). Col. 1: a plume of emission in Pegasus, centred 
at (l,b) = (91?5,-35?8). Col. 2: the Corona Australis region at (l,b) = (0°,-18°). Col. 3: an extended region of emission in 
Musca/Chamaeleon, centred at (/, b) = (305°, -26°) (the Large Magellanic Cloud can be seen to the bottom-right). Col. 4\ the Orion 
(M42) region of emission centred at (/, Z?) = (209°,-19?38). Col. 5: the A Orionis region at (l,b) = (196°,-12°). The graticule 
separation is 10° in both directions, and the colour scales are asinh (although most are close to linear). The regions looked at for the 
correlation analyses are shown as black rectangular boxes. 
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Fig. 14. T-T plots comparing the Commander AME map evaluated at 22.8 GHz to the Commander dust map at 545 GHz in three 
regions, as shown in Fig. 13. The best fit is shown as a dashed line. From left to right, for the Pegasus plume, Musca, and A Orionis 
regions. 


map, which is very closely correlated with the thermal dust am¬ 
plitude at 545 GHz. There is no visible emission in the syn¬ 
chrotron map and thus the signal in the AME map is unlikely 
to be due to synchrotron emission. The feature does have associ¬ 
ated free-free emission, as traced by Her; however, the brightness 
at 22.8 GHz is about 10 times greater than would be predicted 
for free-free emission, assuming a typical electron temperature 
(Te ^ 7000 K). The plume structure is not resolved at 1° resolu¬ 
tion, and therefore absorption of Ha light by dust should be sim¬ 
ilar to other high-latitude regions (\b\ > 30°) at <0.1 mag, with 
only a small fraction being absorbed in compact high-density 
regions. The section of the plume corresponding to S 122 is def¬ 
initely free-free emission; however, we conclude that the rest of 
this feature is predominantly AME, and that its appearance in the 
Commander free-free map is likely to be due to leakage from the 
AME component. This would be an interesting target for higher- 
resolution follow-up measurements, particularly in the bright re¬ 
gion around MBM 53. 

3. Corona Australis. The second column of Fig. 13 shows 
Corona Australis below the Galactic centre (Region 3 in Fig. 10). 
The region has been studied by Harju et al. (1993) in CO emis¬ 
sion and with Spitzer by Peterson et al. (201 1), and it is reviewed 
in Neuhauser & Forbrich (2008). The central object is the R 
Corona Australis dark cloud (R CrA), which has a tail extending 
to the bottom-left of the maps consisting of the refiection nebula 
NGC 6729. The refiection nebula is illuminated by the early-type 
star TY CrA, and exhibits an extended 3.3 p,m emission feature 
that could be to PAH emission (Chen & Graham 1993). We focus 
our analysis on a 3° x 3° area centred on (I, b) = (1°, -22°) that 
encompasses the refiection nebula but not the R CrA cloud. The 
R CrA region is dominated by free-free emission, and it appears 
in both the Ha and free-free maps. However, the refiection neb¬ 
ula clearly shows AME, with a small amount of free-free emis¬ 
sion seen in Ha (there is a feature in the tail in the Commander 
free-free map; this does not show up in Ha, however). The struc¬ 
ture does not appear in the synchrotron map. 

4. p Ophiuchus. This molecular cloud is another 
well-known source of AME (Casassus et al. 2008; 
Planck Collaboration XX 2011). It is included in the PCCS2 
catalogue, so we unmask the point sources when con¬ 
sidering this region. We look at a 5° x 5° patch centred 
on (/,/?) = (353?05,16?9). We find a low AME-lOO^tm 


emissivity of (4.6 + 0.9) p,K(MJy sr“^)“\ compared to (8.3 + 
1.1) pK(MJy sr“^)“^ from Planck Collaboration Int. XV (2014). 
This is due to significant leakage between the AME, free-free, 
and synchrotron components in the Commander solution in this 
region: half of the emission that has been attributed to AME 
in previous analyses (e.g., Planck Collaboration XX 2011) 
is instead attributed to free-free, causing a decrement in the 
synchrotron map, similar to the Perseus region (see Fig. 33 
below). This issue is due to the complexity of this region, which 
lies within the Gould belt, where there is bright surrounding 
free-free and synchrotron emission. 

5. Musca/Chamaeleon region. The third column of Fig. 13 
shows a region of diffuse AME below the Galactic plane, in 
the Musca region close to the Magellanic clouds (Region 5 
in Fig. 10). This region was previously studied in dust polar¬ 
ization by Planck Collaboration Int. XX (2015). Its distance is 
around 160-180 pc, and it contains three dark clouds (Whittet 
et al. 1997). We focus on a 15° x 12° region centred on (/, b) = 
(299°, -16?5). The AME is clearly correlated with the dust emis¬ 
sion, and there is a notable absence of emission in the syn¬ 
chrotron and Ha maps in this region. There is structure evident 
in the Commander free-free map; however, this correlates with 
some parts of the dust emission and not Ha, so it is likely to 
be AME rather than free-free. We mask two bright thermal dust 
sources in the Musca region (Cha I and Cha II) out to a radius of 
80' when calculating the emissivities. 

There is also a much larger AME half-ring that extends south 
of the Musca feature, into the Chamaeleon constellation. We 
include emissivities for the Chamaeleon region in a 35° x 30° 
region centred on (/, Z?) = (305°,-25°). We also use this re¬ 
gion to look at the correlations between AME and IRAS 12 and 
25 p,m, since this part of the sky is at both high Galactic and 
ecliptic latitudes, and as such it is less contaminated by zodia¬ 
cal or stellar light than the rest of the sky. We find an emissivity 
against 12 pm of (550 + 60) pK (MJy sr“^)“^ and against 25 pm 
of (570 + 70) pK (MJy sr“^)“\ both of which are higher than the 
all-sky values given earlier in this section. We also divide the 12 
and 25 pm data by Go (as described above), after which we find 
an emissivity against 12 pm of (250 + 25) pK(MJy sr“^)“\ and 
against 25 pm of (159 + 16) pK(MJy sr“^)“^ That these num¬ 
bers are lower than for the whole sky is likely to be an arte¬ 
fact in the Go map, which has a very small value due to the low 
Commander dust temperature in this region of around 17 K, com¬ 
pared to the mean temperature of 20.9 K. We also correlate the 
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variations in AME emissivities with the fraction of PAH emis¬ 
sion at 12 |Ltm (as described above for the whole sky, following 
Hensley & Draine 2015) in this region. We find similar results to 
above, with a slope of 3.4 + 0.4 when using the Commander ra¬ 
diance map, but poor correlations when using the radiance maps 
based on Planck Collaboration XI (2014). The correlations pre¬ 
sented here are sensitive on how the comparison maps have been 
made, both in terms of inputs to the radiance map and contami¬ 
nating emission in the IRAS maps. 

6. Orion region. The fourth column of Fig. 13 shows the Orion 
region (Region 6 in Fig. 10). The strong free-free emission 
from M42 and Barnard’s Foop (Sh 2-276) is clearly visible. 
However, there is also an arc of AME emission that correlates 
with the high-frequency dust emission. This arc extends to lower 
longitude from M42 along the integral-shaped filament to the 
dark cloud L 1641 and perpendicularly crosses over Barnard’s 
loop (Orion Molecular Cloud A); in the other direction (Orion 
Molecular Cloud B) it extends up to M78 to the top-right of 
the figure. The arc also appears faintly in free-free emission 
and in Ha. This dusty feature was recently mapped in 3D by 
Schlafiy et al. (2015). We look at a 15° x 15° region centred on 
(l,b) = (209?01, -19?38), and we additionally mask M 42 to a 
radius of 100', and the point where the dust crosses the top part 
of Barnard’s loop (where AME has likely leaked into the free- 
free map), when calculating the emissivities. 

7. A Orionis. The A Orionis region (Region 7 in Fig. 10) can also 
be seen to the right of Barnard’s Foop, and is shown in the fifth 
column of Fig. 13. We focus on a 13° x 13° region centred on 
(/, b) = (195?2, -12?2). This Hii region (also known as S 264), 
which is illuminated by the 08 star A Orionis, exhibits a shell 
of AME around the outside of the free-free region coincidental 
with the thermal dust ring. The dust ring has been previously 
noted in Maddalena & Morris (1987) and Zhang et al. (1989), 
and it has been observed in H i by Wade (1957) and CO by Lang 
& Masheder (1998) and Lang et al. (2000). The brightest AME 
regions in the ring are mostly located adjacent to bright ther¬ 
mal dust features, although the brightest AME region does not 
have a bright thermal dust counterpart. AME has been seen in 
PDRs around other Hii regions at higher resolution, e.g., in p 
Ophiuchus (Casassus et al. 2008) and Perseus (Tibbs et al. 2010). 
A Orionis has a particularly high emissivity against 545 GHz and 
T353, but the emissivity against 100 p,m is comparable to the av¬ 
erage; this indicates that the AME is connected to the colder dust 
in this region. It is likely that this region has not been identified 
before due to the large free-free emission feature that the ring of 
AME surrounds. 

There are significant differences between the emissivities for 
these seven regions. Some emissivities (particularly for Perseus 
and p Ophiuchus) are clearly biased low due to component- 
separation issues. Four are relatively consistent with values of 
AME/545GHz typically in the range 40-50 p,K(MJysr“^)“^ 
The notable exceptions are the emissivity values for A Orionis 
and the Chamaeleon region, which are higher than the other 
regions by a factor of 2. However, except for A Orionis and 
Chamaeleon, the emissivities are significantly lower than the av¬ 
erage emissivity of (65 + 7) p,K(MJy sr“^)“^ across the whole 
high-latitude sky. The differences in emissivities could be a 
component-separation artefact, or they could be an indicator of 
the dependence of environmental conditions for the AME, or 


dust grain size distribution. For example, in nearby molecular 
clouds the AME carrier grains could have attached themselves 
to larger dust grains more than in the diffuse medium (e.g., Kim 
et al. 1994). We find a good correlation between the emissivities 
at 545 GHz and the optical depth. 

In conclusion, we find that the Commander AME map pro¬ 
vides a good tracer of AME in our Galaxy, however, there are 
significant degeneracies between the free-free and AME com¬ 
ponents that present difficulties when using the map to calculate 
emissivities. Additional all-sky data at frequencies of 5-20 GHz 
are needed to improve the free-free tracer and so enable a cleaner 
separation of AME from the other low-frequency components, 
in order to determine accurate emissivities and comparison with 
other data sets. In particular, A Orionis is a distinctive region that 
warrants further study. 

4.3. Synchrotron 

In this section we discuss constraints on synchrotron emis¬ 
sion derived from total intensity. We discuss options for mod¬ 
elling the form of the synchrotron spectrum, which is not a 
simple power law, and specify the model used in our baseline 
Commander analysis. The results are strongly limited by un¬ 
certainties in component separation; a much more reliable pic¬ 
ture of the structure of the high-frequency synchrotron emission 
emerges from polarization data, discussed in Sect. 5. However, 
one advantage of total intensity is that we have a high S/N ra¬ 
tio data set in the 408 MHz Haslam map, whereas in polariza¬ 
tion all available ground-based sky maps are at low frequency 
and hence too strongly affected by Faraday rotation to be useful 
for spectral analysis. Because of residual differences between 
Planck and WMAP polarization maps (see Sect. 5.1 and Planck 
Collaboration X 2016), we do not attempt an independent fit of 
the spectrum of the polarized synchrotron emission. 

The Galactic synchrotron spectrum curves significantly be¬ 
low a few GHz (e.g., de Oliveira-Costa et al. 2008; Strong et al. 
2011). To generate a useful spectral constraint from just one 
low-frequency map, we cannot afford to fit this curvature inde¬ 
pendently at each pixel, as otherwise the low-frequency point 
would always fit perfectly and give no constraint on the high- 
frequency spectrum; instead we force the curvature to be con¬ 
stant across the sky, and we have tried several approaches to reg¬ 
ularize the fit. As expected, simple power-law fits are inconsis¬ 
tent with our data, generating large and spurious “gain correc¬ 
tions” at 408 MHz. 

Strong et al. (2011) and Orlando & Strong (2013) model 
the observed synchrotron emission for given Galactic mag¬ 
netic fields and cosmic-ray (CR) propagation scenarios with the 
GALPROP code (Strong et al. 2007). This code solves the CR 
transport equation for any CR species, accounting for diffusion, 
reaccelerating processes, and energy losses in the interstellar 
medium. The CR transport properties are constrained by the lo¬ 
cal CR measurements and the observed ratio of secondary to 
primary nuclei, while energy losses are calculated for a given 
Galactic magnetic field. As implemented in the papers above, 
synchrotron GALPROP models aim to reproduce the large-scale 
emission, but not sub-kpc scale features such as the synchrotron 
loops and individual supernova remnants. The modelled local in¬ 
terstellar energy spectrum of leptons (after propagation effects) 
is adjusted to reproduce the direct measurements by Fermi-LAT 
above 7 GeV (Abdo et al. 2009; Ackermann et al. 2010). At 
lower energies, where solar modulation makes the local spec¬ 
trum hard to determine, the lepton spectrum is based on syn¬ 
chrotron observations. In order to reproduce the curvature of the 


18 


Planck Collaboration: Planck diffuse low-frequency Galactic foregrounds 



0.01 0.10 1.00 10.00 100.00 1000.00 
Frequency / GHz 

Fig. 15. Local spectral index of the synchrotron emission p{v) = 
JlnL/Jlny vs. frequency for a sample of pixels (one per 
A^side = 8 super-pixel), in the GALPROP zl0LMPD_SUNfE model 
from Orlando & Strong (2013). The spectra are colour-coded 
by Galactic latitude: spectra at low latitudes show strong low- 
frequency curvature due to free-free absorption. 


spectrum below a few GHz, Strong et al. (2011) found that the 
injected electron spectrum (before propagation effects) should 
have a break at approximately 4 GeV and should be harder than 
p = 6 below 4GeV. Accounting for all these constraints, the 

resulting injected spectral index used in those works are {N{E) oc 
E~P) with p = 1.6/2.5/2.2, respectively, below/between/above 
the breaks at 4 and 50 GeV. Moreover, the cosmic ray lepton 
(CRL) spectral shape varies with Galactocentric radius due to 
propagation effects, CR source distribution, and the magnetic 
field strength. Despite these effects, the analysis in Strong et al. 
(201 1) and Orlando & Strong (2013) finds that the spectral index 
at any given frequency is strikingly uniform across the sky: the 
spectral index between 408 MHz and 22 GHz has an rms disper¬ 
sion of 0.02; between 22 and 70 GHz the rms is 0.01 (Fig. 15). 

The uniformity of these model spectra over the sky contrasts 
with the much larger variation in p{y) with frequency in the same 
models. For a given CRL spectrum, the frequency of synchrotron 
radiation scales directly with the component of the magnetic 
field perpendicular to the line-of-sight, and therefore the spectral 
index profiles shown in Fig. 15 will be shifted in log frequency 
by the change in \nB^. Given that dfildlnv ^ -0.13 at around 
1 GHz, the spatial uniformity of the model spectral index implies 
that the average value of B^ along each line-of-sight varies by at 
most 15 % (rms) from one line-of-sight to another. Of course, 
each model spectrum is an integral of the emission along the 
line-of-sight, along which Bj^ varies considerably; however, the 
main variation of B is with Galactic scale height z, and hence all 
lines of sight, except those at very low latitude, sample nearly 
the same range of field strengths with nearly the same relative 
weighting. One might expect geometric factors alone to intro¬ 
duce significant spectral variations, but in these models the ran¬ 
dom compo nent of the field dominates, so to a good approxima¬ 
tion ~ V273|5|. 

The synchrotron intensity at fixed frequency depends on B^ 
as jy oc which suggests that regions of enhanced field 

strength will disproportionately dominate the emission. In the 
GALPROP models, the CRTs are supposed to diffuse efficiently 


through the interstellar medium, so that the local CRL density 
is not correlated in detail with the field strength. In reality, if 
the field is enhanced due to compression, e.g., by a shock wave, 
the CRLs will respond adiabatically if the time to diffuse out 
of the region is much longer than the timescale for compres¬ 
sion. Then the particle energies will also be shifted to higher 
energies, giving a larger shift in the piy) profile, roughly cx B^, 
and jy oc ^-2/3-3 exact index depending on the field geom¬ 
etry, see Leahy 1991 for details). Obviously this substantially 
increases the dominance of regions of high B in the overall syn¬ 
chrotron emission. 

The spatial uniformity of the modelled synchrotron spectrum 
is thus likely to be an artefact due to the omission of sub-kpc 
structures in the models, such as the loops and spurs that dom¬ 
inate the high-latitude synchrotron sky. Indeed, we know that 
many individual supernova remnants (SNR) have spectra signif¬ 
icantly flatter than that of the diffuse Galactic emission, as ex¬ 
pected from the steepening caused by energy-dependent diffu¬ 
sion and radiative losses. It may seem redundant to argue on the¬ 
oretical grounds that the GALPROP models show too little spec¬ 
tral structure, since much larger variations have been reported 
based on analysis of existing radio surveys (e.g., Lawson et al. 
1987; Reich & Reich 1988; Dickinson et al. 2009). But if one 
discounts features associated with free-free absorption at very 
low frequencies and free-free emission at higher frequencies, the 
remaining large variations are mostly in the regions of lowest 
synchrotron intensity, which are susceptible to systematic errors 
such as sidelobe contamination from the strong Galactic plane 
emission, ground spillover, zero-level errors, and (at low fre¬ 
quencies) pick-up of sky emission reflected from the ground. 
Reich & Reich (1988) review previous large-area spectral in¬ 
dex maps and conclude that all are unreliable, with the excep¬ 
tion of the 38:408 MHz map of Lawson et al. (1987), which is 
strikingly uniform. Reich & Reich (1988) present their own map 
of spectral index between 408 and 1420 MHz in the northern 
hemisphere showing extensive regions at high Galactic latitude 
with p > -2.4, which they attribute at least in part to free-free 
emission uncorrelated with Her, presumably due to extinction 
of the latter. They later extended their analysis to the full sky 
(Reich et al. 2004), finding even more extensive flat-spectrum 
regions in the southern hemisphere. If this were due to a free- 
free contribution, or even to an unusually flat synchrotron power 
law, the foreground temperature in the steradian or so around 
(/,/?) = (240°,-40°) would be several times higher than ob¬ 
served at 20-30 GHz (e.g.. Fig. la). We believe that the culprit 
is likely to be some combination of the systematic errors men¬ 
tioned above. 

These considerations suggest two ways that our parameter¬ 
ized models can allow for larger spectral variations than found 
by Orlando & Strong (2013). To model emission from regions 
with B significantly different from the average, we can shift the 
spectrum in In v, keeping the overall shape fixed. However, there 
is certainly more spectral freedom allowed than this: from the 
dispersion in the spectral indices of young SNR, we expect the 
injected energy spectrum to be spatially and temporally variable; 
moreover the propagation effects included in the synchrotron 
models of Fig. 15 induce small variations in the high-frequency 
spectrum, which cannot be explained by variations in B. If we 
denote by rGp(v) a fiducial spectrum from that simulation, then 
we can steepen or flatten it by writing 

nyn(v) = (v/vo)"^rGp(y), (4) 

i.e., a change of local spectral index p{y) = JlnT/Jlnv by a 
constant Sp. Most of the variation in Fig. 15 can be accounted 
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Table 3. Key parameters of the WMAP MEM and MCMC fore¬ 
ground models (see Bennett et al. 2013a for details). 


Model 

/^fF 

/^sync 

AME 

MEM 

-2.15 

-3.0 

CNM, Vpeak varies 

MCMC-c base 

-2.16 

-3.0 


MCMC-e sdcnm 

-2.16 

-3.0 

CNM, Vpeak fixed 

MCMC-f fs 

-2.16 

varies 

CNM, Vpeak fixed 

MCMC-g fss 

-2.16 

varies^ 

CNM, Vpeak = 14.95 GHz 


^ As described in Strong et al. (201 1). 


for by such a model, but a larger range of dp is needed to ac¬ 
count for the actual sky. In practice, detailed accounting for the 
spectral curvature below 10 GHz is superfluous given that we 
only use one observed frequency in that regime. For our base¬ 
line Commander analysis our template was the spectrum of a sin¬ 
gle pixel in the zlOLMPD_SUNfE model from Orlando & Strong 
(2013), chosen to be close to the all-sky median for several spec¬ 
tral parameters (i.e., knee frequencies, and spectral indices at the 
knees). As described in Planck Collaboration X (2016), this was 
fitted to the data along with the other foreground components al¬ 
lowing for a global shift in frequency (determined to be a factor 
of 0.26), and an amplitude that was fitted at each pixel. 

The resulting synchrotron amplitude map (shown in Planck 
Collaboration X 2016) is essentially determined by the high S/N 
408 MHz data, and deviates from it only in that residual free-free 
emission at 408 MHz is corrected (or over-corrected, at bright 
compact H ii regions in which free-free absorption is significant 
at 408 MHz). The derived global frequency shift is surprisingly 
large and has the effect of steepening the average spectrum be¬ 
tween 408 MHz and the Planck bands, resulting in a relatively 
low amplitude for synchrotron compared to AME and free-free. 
This may be because the fit is most strongly constrained by the 
North Polar spur and the diffuse halo of the inner Galaxy, the two 
regions that are both strong in synchrotron and relatively free of 
other components that could absorb errors in the synchrotron fit. 
Most analyses find flatter spectra in the narrow Galactic plane 
(e.g., Planck Collaboration Int. XXIII 2015). 

4.4. Comparison with WMAP models 

Here, we discuss the main similarities and differences be¬ 
tween the low-frequency foreground component maps from 
our Commander analysis of the combined Planck and WMAP 
datasets (Planck Collaboration X 2016), and from the final 
WMAP analysis (Bennett et al. 2013a). 

The WMAP team have published sets of component- 
separated maps using a Maximum Entropy Method (MEM) and 
a Markov Chain Monte Carlo (MCMC) technique on a pixel ba¬ 
sis (Bennett et al. 2013b). We have compared the synchrotron, 
free-free and AME maps produced by these methods with the 
Commander maps through T-T plots with HEALPix Aside = 64. 
Table 3 gives the key model parameters of the WMAP models 
that we used. Point sources from the PCCS2 catalogue brighter 
than 1 Jy at 28.4 GHz have been masked out. The best-fit ratios 
of the maps, derived from the T-T plots, are given in Table 4 
for the high-latitude sky, with \b\ > 20° and for \b\ < 20°; a is 
the best-fitting slope. Also listed in Table 4 are the Pearson cor¬ 
relation coefficients (r) for all the comparisons. If r > 0.9, the 
measured slope, a, is reliable, since there is a good linear rela¬ 
tionship between the Commander and the WMAP templates. 


Table 4. Ratios of the WMAP MEM and MCMC maps to the 
Commander maps for the high latitude sky with \b\ > 20° (top) 
and for \b\ < 20° (bottom). A value a > I indicates that the 
Commander maps contain more emission than the WMAP maps. 
The Pearson correlation coefficient (r) is also given. 


Run 

Sync 

Free 

-free 

AME 


a 

r 

a 

r 

a 

r 

MCMC-c base 

0.50 

0.62 

0.68 

0.77 



MCMC-e sdcnm 

0.52 

0.92 

0.77 

0.87 

4.91 

0.75 

MCMC-f fs 

0.52 

0.62 

0.80 

0.77 

3.18 

0.67 

MCMC-g fss 

0.55 

0.62 

0.77 

0.78 

3.14 

0.70 

MEM 

0.34 

0.84 

0.76 

0.79 

2.18 

0.86 

MCMC-c base 

0.12 

0.74 

0.91 

0.83 



MCMC-e sdcnm 

0.61 

0.99 

0.67 

0.98 

2.71 

0.91 

MCMC-f fs 

0.13 

0.78 

1.09 

0.89 

4.71 

0.86 

MCMC-g fss 

0.15 

0.77 

1.06 

0.93 

3.22 

0.88 

MEM 

0.19 

0.90 

1.01 

0.97 

2.36 

0.93 


Slope = 0.5 ±0.1 



0.0 0.2 0.4 0.6 0.8 1.0 

WMAP Synchrotron [mK] 


Fig. 16. Scatter plot between the Commander synchrotron so¬ 
lution evaluated at 22.8 GHz and the WMAP MCMC-f syn¬ 
chrotron model at the same frequency. 


A first obvious difference between the Commander and 
WMAP MEM/MCMC models is that the AME component is 
systematically higher (between 2-4 times) in the Commander so¬ 
lution, at the expense of the synchrotron emission that is lower in 
Commander. This is clearer on the region closer to the Galactic 
plane (bottom half of Table 4), due to the better correlations mea¬ 
sured there. 

When comparing the Commander synchrotron solution with 
the WMAP models that allow the synchrotron spectral index to 
vary over the sky, the data are not well-fitted with a single slope, 
since two populations are clearly present. Fig. 16 shows the scat¬ 
ter plot between the Commander synchrotron solution evaluated 
at 22.8 GHz and the WMAP MCMC-f synchrotron model at 
the same frequency. The two slopes visible in the Figure are 
the result of a flatter spectrum synchrotron component on the 
plane (e.g., Kogut et al. 2007), which is not accounted by the 
Commander synchrotron model, where y^sync is fixed over the 
sky. This is not the case for the WMAP MCMC-e “sdcnm” and 
MEM solutions, which use a fixed synchrotron spectral index. 
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resulting in a much better correlation between the Commander 
and WMAP components (see Pearson correlation coefficients 
for this component in Table 4). We also note from Table 4 that 
the Commander synchrotron component is always lower than the 
WMAP synchrotron models, which is due to the larger AME 
component in the Commander model. 

The Commander free-free component is similar to most of 
the WMAP models at low Galactic latitudes, where Pearson’s 
r > 0.89 for all but the MCMC-c model. The correlation coeffi¬ 
cients in this case are also close to 1. At high Galactic latitudes 
(top half of Table 4), the Commander free-free solution is lower 
than the WMAP models, at the expense of a higher AME con¬ 
tribution in Commander. The wider frequency range of Planck, 
specifically from 100-217 GHz, should enable a cleaner separa¬ 
tion of the free-free component (Planck Collaboration Int. XIV 
2014; Planck Collaboration Int. XXIII 2015). 

In addition to the baseline Commander model, we have also 
compared the WMAP models to several earlier Commander runs 
with different input parameters. In particular, we look at a run 
that used a broken power-law for the synchrotron component 
along with a single AME component, and another that had 
straight power-laws for both the synchrotron and free-free com¬ 
ponents, also with a single AME component. In these models 
the synchrotron spectral index is fitted locally, while the break, 
AyS, in the first is fitted globally. This parameterisation is closer 
to that used by the WMAP team than is our baseline case with 
two AME components and a GALPROP synchrotron spectrum, 
and the fitted results are also closer to the WMAP results, with 
less AME and more synchrotron; however some differences re¬ 
main. In both models the free-free amplitude in the Commander 
outputs is still lower than in the WMAP models, and in both 
Commander still finds less synchrotron emission than WMAP at 
low latitudes (except for WMAP “sdcnm”). However, at high lat¬ 
itudes they find more power in the synchrotron component than 
does the WMAP analysis. The AME components in these ear¬ 
lier Commander models are weaker than in our baseline and in 
much closer agreement with WMAP, with ratios around one for 
both the “fs” and “fss” models at both low and high latitudes, 
while the Commander solutions find more high-latitude and less 
low-latitude emission than the “sdcnm” model. These runs were 
not used for the baseline product as they gave worse overall fits 
and are less physically motivated than our eventual choice; we 
note them here solely to qualitatively consider the effects that 
different models can have on the component outputs. 

The WMAP model that is closest to the baseline Commander 
solution is MCMC-e “sdcnm”, although the AME component 
is brighter in the Commander products. Also, the Commander 
free-free amplitude is about 30% fainter than the correspond¬ 
ing WMAP fit; this has been noted before from comparison with 
RRL data by Alves et al. (2010) and Planck Collaboration Int. 
XXIII (2015). The main limitation in obtaining accurate com¬ 
ponents maps is the precise quantification of the synchrotron 
component. Euskeland et al. (2014) and Vidal et al. (2015) have 
shown that the polarized spectral index of synchrotron emission 
varies across the sky and these variations should be taken into ac¬ 
count to obtain an accurate separation between the synchrotron 
and AME components; low frequency (5-20 GHz) data will be 
crucial for this. 

4 . 5 . LMC and SMC 

The Large and Small Magellanic Clouds (LMC and SMC) 
are satellite galaxies of the Milky Way, close enough 
(50.0 + 1.1) kpc for the LMC, Pietrzynski et al. 2013, and 


(61 + 3) kpc for the SMC, Hilditch et al. 2005) that they are 
resolved by Planck both at 5' in thermal dust and 30' at low 
frequency. They have already been well studied at Planck and 
WMAP frequencies (Bot et al. 2010; Planck Collaboration XVII 
2011; Draine & Hensley 2012). They thus provide a good test¬ 
bed for comparing the Commander solution with previous results 
and expectations, both in the Magellanic clouds and in compar¬ 
ison with our own Galaxy. We show maps of the Commander 
solution in the LMC in Eig. 17, which we discuss component 
by component below. The brightest region in all of the non- 
CMB component maps is the Tarantula nebula, also known as 
30 Doradus, a well-known H ii and star-formation region located 
at (/,/?) = (279?5,-31?7). 

To cross-check the Commander solution, we perform aper¬ 
ture photometry on both the Commander maps and the lat¬ 
est Planck and WMAP maps, as well as other ancillary 
data, to generate spectral energy distributions (SEDs) of the 
LMC and the SMC. Eor the LMC we use central coordi¬ 
nates of (/, Z?) = (279?5,-33?5), with a circular aperture of 
radius 300', and an outer annulus for background subtrac¬ 
tion of 300-350'. Eor the SMC we use central coordinates of 
(/, 5) = (302?8,-44?3), with a circular aperture of radius 150', 
and an outer annulus of 150-200'. We do not attempt to remove 
Galactic foregrounds from the maps or SEDs, except through the 
background annulus in the photometry, since we are using the 
SEDs for the purpose of comparison. As a result, some Galactic 
foreground emission can be seen in the resulting LMC thermal 
dust map in Eig. 17. 

We use the same aperture photometry code that 
was developed for Planck Collaboration XX (2011) and 
Planck Collaboration Int. XV (2014). We have modified the 
code to run on the Commander component maps to measure the 
amplitudes of the individual components in the aperture. The 
uncertainties are derived from the background annulus only: 
they do not include model or calibration uncertainties. Results 
are given in Table 5. We then sum those model SED components 
to compare with the fiux densities measured directly from the 
Haslam, Planck, and WMAP maps, along with radio maps at 
1.4 GHz by Reich et al. (2001) and 2.3 GHz by Jonas et al. 
(1998) and the infrared maps from COBE-DIRBE (Hauser et al. 
1998) and IRAS (Miville-Deschenes & Lagache 2005). This 
comparison for the LMC and SMC is shown in the left panel of 
Eig. 18. We note that the Commander solution does not include 
data points above 857 GHz in its fit. In deriving the models, we 
exclude the 100 and 217 GHz Planck maps, since those bands 
contain CO emission. 

Using the same spectral component model described in 
Planck Collaboration Int. XV (2014), including colour correc¬ 
tions, we also perform least-squares fits (LSE) to the fiux den¬ 
sities up to 3 THz. There are two variants of the LSE model re¬ 
ported in Table 5 and Eig. 18. One is fitted to data including 
the CMB contribution (henceforth referred to as “LSE”) and the 
other is fitted to data post-CMB subtraction, namely fiux densi¬ 
ties from aperture photometry on CMB-subtracted maps (using 
the Commander CMB solution; henceforth referred to as “LSE- 
CMB”). 

We now go through each component in turn for both the 
LMC and SMC, considering both the morphology and the SEDs. 

CMB. The CMB maps for both the LMC and SMC do not ap¬ 
pear to be strongly contaminated by foreground emission; there 
is no clear correlation between the CMB map and the fore¬ 
ground components. The average CMB contribution agrees well 
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Fig. 17. Coimnander solution in the Large Magellanic Cloud region, plotted in Galactic coordinates. Panels show from left to right 
and top to bottom: (a) synchrotron brightness temperature at 30 GHz with 1° resolution (linear colour scale); (b) free-free emis¬ 
sion measure with 1° resolution (logarithmic colour scale); (c) spinning dust brightness temperature at 30 GHz with 1° resolution 
(logarithmic colour scale); (d) thermal dust brightness temperature at 545 GHz with 5' resolution (logarithmic colour scale); (e) 
synchrotron polarization amplitude, P, at 28.4 GHz with 1° resolution (corrected for polarization leakage, linear colour scale); (f) 
CMB temperature with 5' resolution (linear colour scale); (g) CO /=2^1 emission with 5' resolution (linear colour scale); and (h) 
thermal dust polarization amplitude, P, at 353 GHz with 5' resolution (corrected for polarization leakage, linear colour scale). Each 
map covers 15° x 15°, and is centred on Galactic coordinates (/, b) = (279°, -34°). 


in both the Commander and LSF estimates in the LMC and SMC. 
CMB contributions in both objects are positive, as was seen in 
Planck Collaboration XVII (2011). We thus conclude that the 
Commander CMB solution is robust at separating the CMB from 
foreground emission in this region. 


Synchrotron. Low-frequency synchrotron emission is present 
in both the LMC and SMC. It is subdominant in total intensity 
at Planck and WMAP frequencies. In the LMC the peak of the 
emission is slightly offset from the Tarantula nebula, and dif¬ 
fuse emission is also present. The Commander spectral index be¬ 
tween 408 MHz and WMAP/P/a^c^ frequencies is assumed to 
bey^sync = “3.1 according to the GALPROP model (see Sect. 4.3); 
in the LSF and LSF-CMB fits for the LMC we findy^sync ~ -2.7, 
and for the SMC ySsync ~ -2.8, effectively between 408 MHz and 
22.8 GHz. At similar frequencies, Israel et al. (2010) found that 
y^sync = -2.70 + 0.05 in the LMC; however, they found a steeper 
index of -3.09 + 0.10 in the SMC that agrees better with the 
Commander model. As such. Commander will under-estimate the 
synchrotron contribution at Planck frequencies, particularly in 
the LMC and to a lesser extent in the SMC. The synchrotron 
amplitudes at 1 GHz in Commander for both the LMC and SMC 
are significantly lower than that in the LSF, due to the steep 
spectrum assumed, with the difference absorbed by the free-free 
component. The LMC also appears in the synchrotron polariza¬ 
tion map; the polarized emission is at its highest (around 30 % 


polarized, although this would be lower if the Commander syn¬ 
chrotron intensity spectral index were flatter) to the left of the 
peak in the synchrotron total intensity, offset from Tarantula, 
which corresponds to the two polarized synchrotron filaments 
identified by Klein et al. (1993). 


Free-free emission. The free-free component dominates the 
SEDs of both the LMC and the SMC at frequencies of 5- 
50 GHz. The majority of the emission comes from the Tarantula 
nebula; there is also diffuse emission closer to the centre of the 
LMC, and other compact (< 1°) objects are present. The mor¬ 
phology agrees well with the Ha maps of individual sources 
in the LMC by Davies et al. (1976, DEM): the main region 
of 30 Dor is surrounded by a large number of smaller sources 
that show up in the Commander map as diffuse emission. The 
two lower regions are groups of sources, with one compris¬ 
ing DEM 4, 6, and 36 (bottom-left), and the other consisting 
of DEM 27, 28, and 29 (bottom-right). The amplitude of the 
free-free emission can be converted to a star-formation rate 
(SFR) using the equations in Condon (1992). The Commander 
free-free amplitude at 10 GHz is 141 Jy, which gives an SFR 
of O.lOMoyr"^; the LSF model gives 0.08Moyr“\ and the 
LSF-CMB model gives 0.09Moyr“^ These are lower than 
the estimate of 0.2Moyr“^ for the average SFR from anal¬ 
yses of stellar populations by Harris & Zaritsky (2009) and 
Rezaeikh et al. (2014), but agree well with the recent star for- 
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Fig. 18. SEDs of the EMC (top) and SMC (bottom) using aperture photometry, from the Commander solution (left), and least- 
squares fits to the data with (middle, LSE) and without (right)) the CMB (LSE-CMB). The components are synchrotron (orange 
dashed lines), free-free (black double-dashed lines), AME (green dashed lines), CMB (blue dashed lines), and thermal dust (red 
dashed lines). The solid black line represents the sum of the model components. Planck data points are shown in red, WMAP in 
blue, and other ancillary data (some of which are not included in the fits, see text) in black. 
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mation rate of 0.06 M© yr“^ calculated by Whitney et al. (2008) 
based on young stellar objects; these authors also give a range 
of SER estimates of 0.05-0.25 M© yr“^ from infrared and Ha 
data. However, all of these estimates correspond to an SER over 
different timescales (e.g., Murphy et al. 2012), as well as be¬ 
ing subject to systematic effects; a more detailed study would be 
necessary to disentangle these effects. Eor the SMC, the ampli¬ 
tude at 10 GHz of 11 Jy from the Commander and LSE models 
gives an SER of 1.1 x 10“^ M© yr“^; the LSE-CMB model gives 
a slightly higher value of 1.3 x 10“^ M© yr“^ from a flux density 
of 12.8 Jy. These estimates of the SMC free-free amplitude agree 
well with those from Draine & Hensley (2012), who attributed 
11 Jy to the free-free component at 10 GHz, in very good agree¬ 
ment with the Commander amplitude. 


AME. The Commander solution finds a small component of 
AME in both the EMC and SMC. In the EMC, the AME 
map has a bright region centred on 30 Dor, where there are 
many sources seen in the higher-resolution Ha and radio 
maps (McGee & Newton 1972, Davies et al. 1976), and where 
the brightest emission is seen in the EMC in thermal dust. 
The AME towards 30 Dor has a higher peak frequency in the 
Commander solution than that on either side, following the pat¬ 
tern of Hii regions having higher peak frequencies, as noted 
in Sects. 3 and 4.2. At 20 GHz, the peak of the AME com¬ 
ponent, AME contributes about 16Jy to the model, compared 
with the free-free amplitude of around 130 Jy. In the SMC, 
1.6 Jy is attributed to AME, compared with lOJy for the free- 
free component. In both cases, AME is about 10 % of the free- 
free amplitude, and it is comparable to the component sepa¬ 
ration uncertainty in the free-free component, and the ampli¬ 
tudes of the synchrotron and CMB components at that fre¬ 


quency. The ratio of free-free to AME here is significantly higher 
than that seen in the Galaxy (Planck Collaboration XX 2011; 
Planck Collaboration Int. XXIII 2015) due to the presence of 
30 Dor, which dominates the free-free emission in the EMC. 
The LSE analysis includes an AME model, but is consistent 
with little or no AME (the AME contribution seen in the LSE 
SMC SED has a significance below Icr). Nor was AME seen in 
the EMC and SMC SEDs in Planck Collaboration XVII (201 1). 
Draine & Hensley (2012) included an AME component in their 
SED fits to the SMC of 3-5 Jy at 40 GHz; the Commander solu¬ 
tion returns much lower AME amplitudes in the SMC of 0.2 Jy 
at 40 GHz. We calculate AME emissivities for the EMC for a 
15° X 15° region centred on (/, Z?) = (279°,-34°), and for the 
SMC in a 5° X 5° region centred on (/, Z?) = (302?8,-44?3). 
The emissivities are given in Table 2. Eor the EMC we find that 
the AME/545 GHz and AME/T 353 emissivities are comparable to 
those from the AME regions; however, the AME/100 p,m emis- 
sivity is substantially lower, which is likely because the EMC 
has a higher dust temperature than the Galactic average. Eor 
the SMC, we find much lower values for all three emissivities. 
The EMC has comparable PAH levels to our Galaxy (Bernard 
et al. 2008), while the SMC has lower levels, which could ex¬ 
plain why we find higher AME emissivities in the EMC than the 
SMC. The SMC has a very small grain (VSG) population that 
can be seen particularly clearly at 70 p,m (Bernard et al. 2008); 
given the emissivities here this VSG emission is unlikely to be 
connected to AME, as AME is typically thought to be due to 
PAHs. We thus conclude that the small amount of AME found 
by Commander in the EMC is at the level that would be expected 
from our Galaxy, although both the EMC and SMC regions may 
be contaminated by leakage from the free-free and synchrotron 
components, as well as being potentially infiuenced by the steep 
synchrotron spectral index assumed in the Commander model. 
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Table 5. Values of the fitted parameters for the LMC (top) and 
the SMC (bottom) from Commander, and the least-squares fit¬ 
ting before and after CMB subtraction. The Commander param¬ 
eters are as described in Planck Collaboration X (2016) except 
for Async, which has been rescaled to a reference frequency of 
1 GHz. We note that the Commander uncertainties only include 
the standard deviation, and do not include modelling uncertain¬ 
ties. The LSF parameters are the same, except for Ad, which is 
the optical depth at 250 p,m. 


Parameter 

Commander 

LSF (CMB) 

LSF (no CMB) 

EM 

66.2+1.7 

57 + 4 

67 + 2 

Te[K] 

7000 

8000 

8000 

Ad 

24903 + 488 

(3.1 +0.6) X 10-5 

(2.82 + 0.14) X 10-5 

/5d 

1.516 + 0.003 

1.47 + 0.07 

1.42 + 0.03 

Td[K] 

19.29 + 0.09 

22.7 + 0.9 

23.2 + 0.3 

Aamei 

15.3 + 0.3 



AaME2 

3.47 + 0.14 



Acme [qK] 

16 + 3 

15 + 3 


A sync [Jy] 

239 + 4 

406 + 33 

375 + 31 

/^sync 

-3.1 

-2.70 + 0.05 

-2.74 + 0.04 

EM 

20 .0+1.8 

21.7 + 0.9 

25.1+0.4 

Te[K] 

7000 

8000 

8000 

Ad 

2541 + 70 

(1.26 + 0.07) X 10-5 

(7.2 + 0.3) X 10-6 

/5d 

1.463 + 0.004 

1.36 + 0.05 

1.06 + 0.02 

Td[K] 

18.57 + 0.14 

21.5 + 0.3 

25.35 + 0.18 

Aamei 

1.497 + 0.07 



AaME2 

0.57 + 0.03 



Acme [qK] 

36 + 5 

37 + 2 


A sync [Jy] 

25 + 2 

34 + 3 

31 + 3 

/^sync 

-3.1 

-2.76 + 0.11 

-2.88 + 0.11 


Further work is needed to improve constraints on the AME and 
its properties in this region. 

CO. The Commander CO /=1^0 map detects various regions 
of CO emission, particularly around the 30 Dor region, as well as 
in various other dusty regions in the LMC. The morphology of 
the map compares well with ground-based surveys, such as the 
Magellanic Mopra Assessment (Wong et al. 2011). This indi¬ 
cates that the separation of CO emission has worked well in this 
region. An exception is a large negative region in the Commander 
high-resolution map, which is a component separation artefact. 

Thermal dust. The thermal dust seen by Planck agrees well 
in morphology with IRAS data. The dust temperatures from 
Commander are significantly lower than those from the LSF: 
for the LMC the average dust temperature in the aperture is 
(19.29 + 0.01) K, compared with (22.7 + 0.9) K for the LSF. 
Planck Collaboration XX (2011) found (21.0 + 1.9) K when fit¬ 
ting for a free (3^ = 1.48 + 0.25 (comparable with the (3^ found 
here). Commander only fits data up to 857 GHz and the fitted 
temperature and spectral index underestimates the flux densities 
at higher frequencies. The LSF fits data up to 3 THz and so finds 
a higher temperature. We also note that where the fitted tem¬ 
perature is higher, (3^ is lower. The same applies for the SMC, 
where the average from Commander is (18.57 + 0.14) K, com¬ 
pared to (21.5 + 0.3) K from the LSF, and (22.3 + 2.3) K from 
Planck Collaboration XX (2011), where (3^ = 1.21 + 0.27 is fiat- 
ter than found by Commander but in better agreement with the 


LSF.^ In the SMC, the submm excess (see e.g., Israel et al. 2010, 
Planck Collaboration XVII 2011) has been subsumed into the 
dust spectral index and the lower dust temperature. The thermal 
dust polarization in the LMC traces the spiral structure, with the 
projected magnetic field running parallel to the structure, partic¬ 
ularly in the 30 Dor region where the Planck data have a higher 
sensitivity due to its scan strategy. This is also where the syn¬ 
chrotron polarization is seen. 

In conclusion, we find that the Commander component maps in 
the regions of the LMC and SMC largely agree well with previ¬ 
ous results and expectations, although the thermal dust proper¬ 
ties are only representative up to 857 GHz. 

5. Polarized foregrounds at 1 cm 

The previous section showed that significant uncertainties re¬ 
main for component separation in total intensity. Fortunately, the 
picture is quite different in polarization, where, between Planck 
and WMAP, we have twelve bands with maps of the all-sky po¬ 
larization. Our current understanding is that there are three sig¬ 
nificantly polarized components: the CMB; synchrotron emis¬ 
sion; and thermal dust. AME polarization has not been clearly 
detected, with upper limits of typically a few per cent (see 
Sect. 5.7). 

The three strongly polarized components have radically dif¬ 
ferent spectra, and moreover, CMB polarization is well separated 
in angular scale from the foregrounds: after convolution to 1° 
resolution, the CMB polarization predicted by the Planck best- 
fit cosmology contributes an rms of only 0.54 yuK to the Q and U 
maps, negligible compared to the noise per beam (although, of 
course, detectable in the angular power spectrum at low multi¬ 
poles). The spectrum of the polarized dust is discussed in Planck 
Collaboration Int. XXII (2015) and Planck Collaboration Int. 
XXX (2016). Here we are interested primarily in the synchrotron 
component and consider the dust polarization only to the extent 
that the magnetic field pattern that it traces sheds light on the 
synchrotron features (Sect. 5.2). Planck Collaboration Int. XIX 
(2015) informs our understanding of the dust polarization fea¬ 
tures. 

We begin by constructing a new map of polarized syn¬ 
chrotron emission by combining WMAP and Planck maps, 
which significantly increases the average S/N ratio compared 
with Planck data alone. We then discuss some of the 

major features in the polarized sky, including the loops, spurs, 
filaments, and bubbles. 

5.1. Combination of Planck and WMAP 

WMAP (Page et al. 2007; Bennett et al. 2013a) pro¬ 
vided the first clear view of the intrinsic synchrotron po¬ 
larization across the sky - previous ground-based observa¬ 
tions (e.g., Brouw & Spoelstra 1976; Wolleben et al. 2006) be¬ 
ing strongly affected by Faraday rotation and depolarization 
due to their lower observing frequencies (< 2 GHz). WMAP’s 
lowest-frequency channel, K-band, dominates the combined 
WMAP S/N ratio in polarization because of the steep spectrum 

^ The Commander solution has a Gaussian prior on the thermal dust 
spectral index of = 1.55 + 0.1, and the fitted values for the LMC 
and SMC are towards the lower end of this prior; however Commander 
would have returned a lower value for (3 than the prior if the data had 
preferred it. The LSF indicate flatter spectral indices, however this is 
including data over a larger (> 857 GHz) range of frequencies. 
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Fig. 19. Difference between the combined weighted polariza¬ 
tion maps for Planck and WMAP at 1° resolution, defined as 


[(Gpianck “ GwMAp)^ + (f^Planck “ f^WMAp)^] 


1/2 


of synchrotron emission. Although the full-mission Planck-LFl 
maps have significantly lower noise than the final WMAP re¬ 
sults, the foreground S/N ratio is very similar, since at Planck's 
lowest frequency (28.4 GHz) the synchrotron brightness is half 
that in K-band. After smoothing to a common resolution of 
1°, the median (mean) S/N ratio is 2.47 (3.77) for WMAP K- 
band and 2.64 (3.72) for Planck 28.4 GHz. The different scan 
strategies of the two missions result in somewhat different sen¬ 
sitivity patterns, so that each map is superior in some sky re¬ 
gions; Planck has the largest advantage in the regions around the 
Ecliptic poles (see the hit count maps in Planck Collaboration VI 
2016). 

Planck and WMAP polarization results generally agree well, 
but the difference map (see Fig. 19) shows significant artefacts, 
both on the plane and at high latitudes. Extremely large angular- 
scale {I = 3,5,7) residuals covering much of the sky are thought 
to be due to the poorly-constrained modes in the WMAP data 
as a result of the scan strategy (Jarosik et al. 2011; Bennett et al. 
2013a); they occur in all the WMAP frequency maps, with an ap¬ 
proximately consistent pattern. However, they are most readily 
visible in the V- and W-bands of WMAP when the foreground 
brightness is close to the minimum. The on-plane features are 
likely to be caused by inaccuracies in the Planck leakage maps. 
These differences are further discussed in Planck Collaboration 
X (2016). 

Because polarization is only marginally detected in many 
pixels of the individual WMAP and Planck maps, we have con¬ 
structed a weighted mean image using all the maps in the range 
20-50 GHz, where the polarization is dominated by optically 
thin synchrotron radiation, with negligible Faraday rotation (ex¬ 
cept at the Galactic centre, Vidal et al. 2015), so that the polar¬ 
ization angles should be consistent between the maps. We first 
smooth the Planck and WMAP maps to a resolution of 1° and 
Aside = 256 in both Q and f/, and the QQ, UU, and QU covari¬ 
ances as well. We assume that the emission has a single spectral 
index of -3.0, and we use this to scale all the maps (and their 
covariances) to match the Planck 28.4 GHz image. We then use 
the 2 X 2 2 and U covariance matrices and their inverse ma¬ 
trix to calculate the weighted average in each pixel of the Q 
and U maps. We create three weighted maps as a result: one 
consisting of the three lowest WMAP bands; one consisting of 
the two lowest Planck bands; and one using both WMAP and 


Planck. The Planck variance estimates include the uncertainty 
in the bandpass leakage correction. This results in a large frac¬ 
tional error in Q and U along the plane, so WMAP data dom¬ 
inate the on-plane emission in the combined map. At high lati¬ 
tudes, the somewhat higher sensitivity of Planck down-weights 
the impact of the large-scale WMAP polarization artefacts, but 
we note in our discussion when these may still have an impact 
on our results. The improvement in S/N ratio is significant, with: 
a median (mean) S/N ratio of the WMAP weighted map of 2.70 
(4.19); the Planck weighted map of 2.68 (3.81); and the com¬ 
bination of both of 3.37 (5.39), where the bigger improvement 
comes from the combination of both WMAP and Planck. This 
corresponds to a roughly 25 % improvement in S/N ratio com¬ 
pared to using the combined WMAP or Planck maps separately. 

Commander also provides a polarized intensity synchrotron 
map, which is calculated using only Planck data. As this map is 
the best-fit solution to a model using noisy data, it is noisier than 
our combined Planck-on\y map (our map treats the small CMB 
component as noise). Moreover, the addition of WMAP K-band 
data in our map further reduces the noise in the final combined 
map. At 1° resolution, the Commander map has a median noise 
value over the entire sky of (Tq ^ crjj = 4.5 jjK, while our com¬ 
bined map has ctq ^ ctu = 2.6yuK. 

Polarized intensity maps have been debiased using the 
asymptotic estimator (Montier et al. 2015; Vidal et al. 2014), 
which generalizes the estimator first proposed by Wardle & 
Kronberg (1974) to the case of anisotropic errors in (Q, U). 
Figure 20 shows the resulting polarization combination maps. 
A number of new polarized structures visible in the combined 
map that were unclear in the individual maps are highlighted. 
These are discussed in the sections that follow. 

5.2. Overview of polarized synchrotron emission 

Synchrotron total intensity is distributed comparatively uni¬ 
formly over the sky. In the 408 MHz map, assuming an instru¬ 
mental plus extragalactic background of 8.9 K (Wehus et al. 
2014), half the total Galactic flux is contributed by 18 % of the 
sky.^^ The equivalent figure for the dust-dominated Planck 545- 
GHz map is 4.6 %. In the synchrotron intensity map, although 
the Galactic plane and the North Polar spur are visually promi¬ 
nent, they are superposed on broad diffuse emission that domi¬ 
nates the total flux. 

This diffuse emission seems to be much weaker in polar¬ 
ization. Away from the narrow Galactic plane, i.e., at \h\ > 3°, 
the polarized cm-wavelength emission is dominated by the syn¬ 
chrotron loops and spurs familiar from low-frequency radio sur¬ 
veys. In fact, Vidal et al. (2015) demonstrate a close correspon¬ 
dence between the polarized intensity at WMAP K-band and an 
unsharp-masked version of the 408 MHz map, in which structure 
on scales > 10° is Altered out (see Fig. 21 top row a and b). 

In particular, there is hardly any trace in the Planck and 
WMAP polarized maps of the synchrotron halo of the inner 
Galaxy, which Alls roughly |/| < 60° and 5° ^ \b\ ^ 15° 
(Planck Collaboration Int. XXIII 2015) and contributes 20-30 % 
of the Galactic synchrotron flux in low-frequency maps (com¬ 
pare panels a and c in Fig. 21). The effect is even clearer in 
fractional polarization, which is very low (< 10 %) in the in¬ 
ner halo region (see Fig. 22). The component separation anal- 


The diffuse emission may be even brighter, since the background 
assumed may include an isotropic component of Galactic emission. 
Using the Lawson et al. (1987) extragalactic background estimate of 
5.9 K, half the sky flux comes from 21 % of the sky. 
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Fig. 20. Top: Combined weighted polarization intensity map after debiasing, with features highlighted. The black dash-dot lines 
show the outlines of Loops I to IV, as defined by Berkhuijsen et al. (1971), the blue dashed lines show the filaments described by 
Vidal et al. (2015) using WMAP polarization data, the red dashed lines show features that are visible in the new Planck data, and the 
magenta dashed lines show the outline of the Fermi bubbles. Filaments that are discussed in the text are labelled. Bottom: The same 
combined polarization intensity map, with projected magnetic field angle (at 90° to the polarization angle) encoded in colour with 
asinh scaling. The coloured half-disc represents the polarization angle depicted in the map, with 0° horizontal, while the polarization 
intensity is represented with the radial distance along the half-disc. 


ysis of Sect. 4 implies that the halo is still present in the un¬ 
polarized cm-wavelength maps, because the data are well-fitted 
with our synchrotron model, which has a constant shape, so that 
the synchrotron distribution is no different at 20 or 30 GHz from 
408 MHz. Although, as we have seen, this separation is subject 
to substantial uncertainties, it is corroborated by spectral analy¬ 
sis of ground-based surveys (Reich & Reich 1988; Platania et al. 
2003) that find that the halo and spurs have similar spectral in¬ 
dices near 1 GHz. 


The magnetic field in the Galactic disc follows the plane, as 
expected from the shearing effect of differential rotation, and one 
would expect the same to apply in the inner halo. Its polarization 
would then be roughly orthogonal to that in the spurs rising out 
of the plane, which all have fields roughly parallel to their axes 
(Vidal et al. 2015). Therefore, the net polarization of halo and 
spurs will cancel to some extent. But the fact that the spurs are 
prominent in polarization, even when they are only perturbations 
on the inner halo in total intensity (see Fig 23), implies that the 
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Fig. 21. Various tracers of the interstellar medium in the hemisphere containing Loop I, with all maps centred at I = 330°, b = 0°, 
and using stereographic projection, so that small circles in the sky are projected as circles: (a) 408 MHz; (b) unsharp-mask 408 MHz; 
(c) combined Planck-WMAF polarization intensity; (d) Hi at vlsr = 0, velocity width 10kms“\ from the LAB survey (Kalberla 
et al. 2005); (e) the same, at -10kms“^; (f) Planck 353 GHz temperature; (g) Planck 353 GHz polarization intensity; (h) ROSAT 
bands R4-rR5 (0.44-1.21 keV) from Snowden et al. (1997); and (i) Fermi 1-10 GeV from Ackermann et al. (2014). The dashed 
white outline is our proposed outer boundary, derived from panels (b) and (c). 
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0.0 ■ 70.0 % 


Fig. 22. Polarization percentage of the debiased weighted polar¬ 
ization intensity map against the Commander synchrotron total 
intensity fit evaluated at 22.8 GHz. Uncertainties in the compo¬ 
nent separation discussed in the text, and also uncertainty in the 
absolute zero level, make this highly uncertain, except in regions 
of strong synchrotron emission, such as the North Polar spur and 
the inner Galactic halo (see also discussion in Vidal et al. 2015). 


halo must be very weakly polarized; even between the spurs, 
there is hardly any sign of a field parallel to the plane in the in¬ 
ner halo region (a trace may be visible near (/, b) = (34°, 8°); see 
the discussion of the I = 45° feature in Sect. 5.4). This is a sur¬ 
prising contrast to the narrow plane at |Z?| < 3°, where the overall 
parallel orientation of the field is very clear in Fig. 20, despite 
the fact that a significant fraction of its emission is due to in¬ 
dividual SNRs (Planck Collaboration Int. XXIII 2015), whose 
overall polarization orientation is nearly random. The conven¬ 
tional analogy between the Galactic halo and the Solar corona 
suggests that the halo field should relax to a nearly force-free 
configuration with little small-scale structure, constrained pri¬ 
marily by the foot points where the field lines are tied to dense 
gas clouds in the plane. This would lead us to expect a reason¬ 
ably high fractional polarization; differential rotation of the foot- 
points should, as with the disc field, shear the field so that it is 
largely parallel to the plane as viewed from Earth. Evidently the 
halo field is much more tangled than this naive argument would 
suggest. 

It is worth noting that the overall fractional polarization also 
appears to be low (< 15 %) across the high latitude sky (Fig. 22), 
although this quantity is very sensitive to the poorly-known zero 
level of the synchrotron total intensity (Vidal et al. 2015). Again, 
the implication is that there is substantial tangling of the field 
even on lines of sight looking out of the disc. 

5.3. Loop I 
5.3.1. Structure 

Loop I is the nearly-circular structure of radius 58° whose top- 
left quadrant (as viewed in Galactic coordinates)^^ is traced 
by the North Polar spur (NPS). It is also detected in soft X- 
rays, where the emission is dominated by thermal emission from 


To avoid the non-intuitive concepts of Galactic “East” and “West”, 
throughout this section we describe emission features in terms of left 
and right as projected in the figures, corresponding to the directions of 
increasing and decreasing longitude, respectively. 



Rodius [°] 


Fig. 23. Profiles of Loop I in total intensity at 408 MHz 
(Remazeilles et al. 2015) and in our combined Planck-WMAF 
polarized intensity image at 28.4 GHz. Profiles are along ra¬ 
dial cuts beginning at the nominal centre of the loop, (/, b) = 
(329°, 17?5) and are averaged over a 10°-wide sector in position 
angle, at position angles 75° (solid line), 50° (dotted), and 25° 
(dashed). They cross the North Polar spur at Z? = 22°, 43°, and 
61°, respectively. 


3x 10^ K gas (Willingale et al. 2003), and is bordered by cold ma¬ 
terial visible via H i and dust emission (both thermal and anoma¬ 
lous). The NPS has long been a suspected y-ray emitter, and 
Ackermann et al. (2014) clearly detect Loop I at GeV energies, 
presumably due to inverse-Compton scattering of starlight by the 
CRLs in the loop, combined with pion decay emission from the 
cold border. Figure 21 shows the structure in several tracers in 
stereographic projection (chosen because circles are projected 
as circles and also because the angular scale increases outwards, 
enabling a clear display of features around the edge of this enor¬ 
mous structure that covers most of a hemisphere). 

In both Galactic hemispheres the polarization maps show a 
number of spurs within Loop I that parallel its outer boundary, 
nearly all concentrated on the left-hand side of the structure. 
Since Loop I covers about a third of the sky, including the in¬ 
ner Galaxy, some of the features projected inside it are surely 
unrelated, but the general coherence of the structure strongly 
suggests that most of the emission away from the plane has a 
common cause. The internal spurs are much more obvious in 
polarization than in total intensity, even after unsharp masking 
(Figs. 21a,b,c, 23). These features are therefore highly polar¬ 
ized; however, they are projected onto diffuse extended emission 
that seems to be weakly polarized (< 15 %), including the inner 
halo of the Galaxy, so the overall fractional polarization towards 
the inner spurs is not much higher than towards the NPS itself 
(~ 40-50 %), while between them the fractional polarization is 
low (Fig. 22). 

Previous measurements of the radio outline found Loop I 
to be surprisingly close to circular: the definitive study by 
Berkhuijsen et al. (1971, hereafter BHS) found that 19 points 
on the ridge-line covering 155° around the loop (all in the north 
Galactic hemisphere) fitted a small circle with an rms of only 
0?9. The NPS conforms roughly to the brightness profile ex¬ 
pected from a shell of emission, for which the ridge-line marks 
the tangent to the inner surface; the outer boundary lies several 
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degrees beyond the ridge line, and is also clearly traceable, espe¬ 
cially around the NPS (e.g.. Fig. 21b). Relatively bright diffuse 
polarized emission outside the boundary of the NPS south of 
b = 40° (Fig. 21c, 25) seems to be part of a different structure 
(see Sect. 5.4). 

South of the plane, the loop is superimposed on the inner 
halo and cannot be clearly followed in total intensity, but in 
WMAP and Planck polarization maps (Figs. 20,21), the pat¬ 
tern is clear, at least on the left-hand side of the structure. The 
brightest southern spur, called filament Is by Vidal et al. (2015), 
runs well inside the path of the BHS small circle. However, a 
fainter spur extends from (/,/?) = (22°,-5°), closely along the 
BHS circle until its ridge becomes lost in faint diffuse emis¬ 
sion near b = -37°; at this point clearly-detectable polariza¬ 
tion again extends several degrees outside the circle. There is 
then no trace of a ridge-line close to the BHS circle along its 
whole bottom-right quadrant. It eventually re-aligns with the 
observed loop boundary at the top of the spur north of the 
Vela SNR (/,/?) = (268°, 25°), near the first measured point 
in Berkhuijsen et al. (1971). In this bottom-right quadrant three 
broad, faint, concentric arcs extend from the interior of Loop I, 
reaching some 30° beyond the BHS circle; these arcs are vis¬ 
ible both in polarization and in the unsharp-masked 408 MHz 
map. The outermost (filament XII of Vidal et al. 2015) rejoins 
the plane near I = 250° and appears to emerge on the other side 
as the aforementioned spur above Vela. Although they might be 
unrelated, they share the general pattern of the other Loop I fil¬ 
aments, for instance the middle arc that passes in front of the 
LMC is clearly brightest on the left-hand side, as it approaches 
the Galactic plane near I = 300°. We consider it plausible that 
filament XII represents the outer boundary of Loop I, and on 
this basis we have marked the outline of the loop in Fig. 21. 
None of the current theories (see Sect. 5.3.3) predict a strictly 
spherical structure; taking the smallest plausible distance to the 
centre as 120 pc, the diameter is at least twice the 100 pc scale 
height of the H i distribution. Therefore, even ignoring the pos¬ 
sible interaction with the Local Cavity (LC) in which the Sun 
is situated, the loop is expanding into an inhomogeneous envi¬ 
ronment. Most individual SNRs, even though much smaller than 
Loop I, show quite large departures from circular outlines, often 
exceeding the level of distortion implied if filament XII traces 
the edge of Loop 1. 

In X-rays, the North Polar Spur peaks at a position inwards 
from the synchrotron ridge, but, as the outline drawn in Fig. 21 
shows, the outer edge of the spur in the two tracers is coincident. 
This result is unexpected, because the Spur cannot be bounded 
by a shock front fast enough (v > 300 km s“^) to heat the ambient 
medium to millions of kelvin, while at the same time accelerat¬ 
ing the cold border to at most 25 km s“^ (see Sect. 5.3.3 below). 
In young supernova remnants, co-incidence of radio and X-ray 
boundaries is always due to X-ray synchrotron emission from 
the shock, and it seems likely that in the NPS the faint X-ray 
emission from the outer edge is also non-thermal. 

In contrast, the gas and dust tracers extend a few degrees 
further out, particularly obvious in the H i and dust filament at the 
top of the loop, (/, b) = (320°, 84°). This is just as expected, since 
neutral atoms and dust grains could not survive at the 3 x 10^ K 
of the X-ray emitting gas. 

5.3.2. Distance 

Loop I is usually associated with the Sco-Cen OB association 
(e.g., Salter 1983) at a distance D = 120-140 pc (de Zeeuw et al. 
1999). The primary evidence is that the NPS appears to be de¬ 


tectable in starlight polarization for stars at distances of > 100 pc 
(Mathewson & Ford 1970; Santos et al. 201 1); the aligned grains 
must be in the cold border rather than the spur itself. Planck-Wl 
maps of polarized dust (Fig. 21g and Planck Collaboration Int. 
XIX 2015) show that emission and extinction measurements of 
field direction agree, confirming that the extinction distances ap¬ 
ply to the dust features seen in emission. 

Sofue (1977,1994) has instead suggested that the loop is due 
to an outburst at the Galactic centre, but in Sect. 5.5 we argue 
that this model is now ruled out. However, evidence against the 
traditional distance has also been accumulating. 

I wan (1980) showed that the wall of dense gas that bounds 
the LC would have to form a major “dent” in the apparently 
spherical structure of Loop I, assuming D = 130 pc. Placing the 
loop several hundred parsecs away avoids the awkwardness that 
this dent is not reflected at all in the projected outline. If the 
shape is quasi-spherical, its near surface would only be about 
20 % of the distance to the centre, D ^ 400 pc given absorp¬ 
tion at the LC wall. It would still be plausible that stars at about 
100 pc could be polarized by dust in the border. 

The interstellar medium within 200-300 pc has now been to- 
mographically mapped in a variety of absorption tracers against 
stars with Hipparcos distances (Frisch et al. 2011). Towards the 
bright left side of Loop I, the LC wall is at D = 70-80 pc, 
and in the Galactic plane dense gas then extends continuously 
to the Ophiuchus and Lupus molecular clouds at D = 140- 
160 pc (e.g., Vergely et al. 2010). These are the remnants of the 
parent cloud(s) of the Sco-Cen association. Although there is 
a small cavity surrounding the most active part of the associa¬ 
tion (Puspitarini & Lallement 2012), there is no sign of a cav¬ 
ity with the angular size of Loop I, and this can no longer be 
ignored given that microwave polarization maps show that the 
loop crosses the Galactic plane without interruption. Lallement 
et al. (2014) push the analysis out to 1 kpc by using stars with 
photometric distances instead of relying purely on Hipparcos. 
They find a low-density cavity starting some 250 pc from the 
Sun, covering the longitude range of Loop 1. If this is the Loop, 
the radius is roughly a kiloparsec, and it would stretch well into 
the halo. 

Given the lack of room for Loop I immediately 
outside the Local Cavity, Frisch has long argued 
(e.g., Frisch 1981; Frisch et al. 2011) that the two are con¬ 
tinuous, and that the NPS is merely separated from us by a fold 
in the LC wall. However, this view is at odds with the evidence 
from soft X-ray absorption of the NPS. This absorption is 
visible in the ROSAT map (Fig. 21h) where the inner boundary 
of the X-ray NPS at Z? < 40° closely follows the border of the 
cold material traced in panels d-g of Fig. 21. This material can 
be confidently assigned to the LC wall, since it is picked up 
by starlight polarization at D > 100 pc (Santos et al. 2011), 
and again extinction and emission polarization is consistent. 
X-ray spectra in the region are well fitted by absorbed thermal 
emission (Puspitarini et al. 2014). Similarly, Sofue (2015) 
makes a strong case that the sharp truncation of the X-ray NPS 
below b = 10° is due to absorption in the cold clouds of the 
Aquila Rift. Both papers find that the hydrogen column density 
derived from the X-ray spectra requires a distance to the front 
of the spur, roughly the shell tangent-point, of at least 200 pc. 
The corresponding distance to the centre is D > 400 pc. 


Sofue (2015) argues for a distance of ^ 500 pc, but this relies on the 
assumption that the Aquila Rift is at a single distance, whereas clouds 
at several different distances are involved, see the tomographic maps of 
Lallement et al. (2014). 
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Fig. 24. Line elements showing the orientation of the projected magnetic field (90° to the synchrotron polarization angle) from the 
Planck 28.4 GHz data at 2° resolution. The region within b = +10° has been masked out (grey rings) to avoid crowding of the 
vectors. The colour scale shows the polarized intensity. The maps are in orthographic projection, centred on the Galactic poles, 
North (left) and South (right). Light dashed circles show the outlines of Loops I to IV, as defined by Berkhuijsen et al. (1971) and 
shown in Fig. 20; the red dashed line is the locus of filament Ills from Vidal et al. (2015). 


Further evidence for a larger distance comes from the 21- 
cm polarization map of Wolleben et al. (2006). As noted by 
Wolleben (2007), the emission from Loop I is strongly depo¬ 
larized at \b\ < 30°. Given the cut-off in latitude, which does not 
correspond to any feature in the intrinsic intensity or polarization 
structure, the depolarization is almost certainly caused by fluctu¬ 
ations in the foreground Faraday depth, (pp, across the 30-arcmin 
beam of the Wolleben map. The sharp cutoff in b suggests that 
the high-latitude part of the loop projects above the dense layer 
of the ISM that creates the “Faraday horizon”. At the near dis¬ 
tance, the implied scale height, h ^ 50 pc, is much too low for the 
most effective source of widespread Faraday rotation, the warm 
ionized medium (WIM) {h = 1-2 kpc, Gaensler et al. 2008). 

The Faraday depths required also imply a larger distance. 
Depolarization requires fluctuations with crp ^ = 

36radm“^. The LC has negligible Faraday depth, due to its 
low density, and so the magneto-ionic medium responsible must 
be in the LC wall, if Loop I is directly adjacent to us, even 
though the wall is mainly traced by neutral gas. Fluctuations 
across a beamwidth 6 requires field tangling on scales DO, and 
hence the Faraday depth along the line of sight of length L cen¬ 
tred at distance D integrates as a random walk, with amplitude 
<^RM = ^F(^el^iil) ^/LDO, where kp = 8.1 x 10“^ mnT“^ pc“^ 
The WIM has electron density ^ 1 cm~^ at pressure equilib¬ 
rium. We also assume B = I nT, with (5||) a factor 1/ V3 smaller. 
Inserting these values and taking D = 80 pc, the required path 
length is L ^ 70 pc, which places the front face of the loop be¬ 
yond the the Sco-Cen association, even if the WIM has unity 
filling factor. Of course higher-density material would require 
less path and in fact the Sh 2-27 H ii region (G006+236), ionized 
by f Oph at D = 112 pc (van Leeuwen 2007), casts a partic¬ 


ularly deep Faraday depolarization “shadow” at 21 cm, imply¬ 
ing it is in the foreground despite a precise distance that puts it 
near the assumed shell centre. However, such objects are readily 
detected in Ha surveys, and cover only a small fraction of the 
depolarized region. The near distance can only be made com¬ 
patible with these results if the polarization is dominated by the 
far hemisphere of the loop, as in one of the models discussed in 
Sect. 5.3.3. 

Given all this evidence that Loop I is several times further 
away than the Sco-Cen OB association, we searched for an al¬ 
ternative group of OB starsto generate the Loop. Given the 
angle at which the NFS passes through the Galactic plane, these 
must be at positive latitudes; we searched -45° < I < 15°, 
5° < b < 35°. Apart from Sco-Cen, in this region there are 
no known OB associations (de Zeeuw et al. 1999), but an asso¬ 
ciation old enough to generate Loop I could already have lost its 
O stars and could easily have escaped notice. Using SIMBAD 
(Wenger et al. 2000) we found 51 OB stars (mostly B giants) 
with parallax, tt, between 1 and 3.3 mas; 14 with 0 < tt < 1 mas; 
and 14 with tt < 0. Although the individual parallaxes are barely 
significant, the distribution is strongly biassed to positive values, 
so there certainly has been star formation within the last 20 Myr 
in the relevant volume. We will have to await Gaia data for a 
definitive assessment of stellar groupings there. 


We follow Reed (2003) in defining OB stars as main sequence stars 
of type B2 and earlier, and giants of type B9 or earlier, roughly corre¬ 
sponding to stars massive enough to go supernova. 
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Fig. 25. Combined 28.4 GHz polarization map at 1° resolution of 
the North Polar spur and the region outside its perimeter, includ¬ 
ing filament X of Vidal et al. (2015), which leaves the Galactic 
plane at a shallow angle from near I = 36°, and the I = 45° fea¬ 
ture between filament X and the NPS. Colour scale is asinh, but 
the length scale for polarization line elements remains linear, as 
in all our figures. 


5.3.3. Interpretation 

The majority view of Loop I is that it is a nearby pre-existing 
cavity re-energized by one or more recent supernova, in order to 
explain the cold neutral border seen in dust and H i (see Salter 
1983 for a review). 

The Hi border appears at Vlsr = 0kms“^ (Fig. 2 Id), which 
is uninformative, because expansion at the tangent points is per¬ 
pendicular to the line of sight. Heiles (1984) interpreted the 
border as part of a shell expanding at 25kms“^; this im¬ 
plies that the obvious Sco-Cen supershell, visible most clearly 
at -10kms“^ (e.g., Kalberla et al. 2005; Vidal et al. 2015, 
Fig. 21e), is merely the end-cap of Loop I itself. It is worth not¬ 
ing that this shell is elongated in the same direction as our pro¬ 
posed outline for Loop I. However, the Sco-Cen shell is more 
usually seen as a distinct structure with a smaller angular size 
than Loop I (e.g., de Geus 1992). If Loop I has a much lower 
expansion speed, e.g., 2kms“^ as estimated by Weaver (1979), 
its HI emission would not be separable in velocity from the very 
local HI surrounding the Sun. The difference between these two 
interpretations is dynamically significant: at 25kms“^ the cold 
shell cannot have been accelerated by a shock, since the heating 
would have dissociated the gas and evaporated the dust; there¬ 
fore acceleration by a pressure gradient as in the stellar wind 
model proposed by Weaver (1979) is required. A 2kms“^ ex¬ 
pansion is consistent with the weak remains of a shock, although 
of course a much faster shock would be needed to heat the inte¬ 
rior gas above 10^ K. We also note the absence of optical emis¬ 
sion lines characteristic of a cooling shock, which would be ex¬ 
pected if the shock speed exceeded 10kms“^ A plausible ex¬ 


planation is that the re-energizing blast wave has only recently 
hit the cavity wall, which would also explain why the NPS is 
so much brighter and more sharply-defined than the other large 
loops (e.g., Borken & Iwan 1977). 

Figure 24 shows the projected magnetic field orientation 
around the Galactic polar caps. We show the Planck data only, 
smoothed to 2° resolution, because of the suspect large-scale 
structure in the WMAP data. Figure 25 shows the pattern closer 
to the Galactic plane. The quasi-parallel field pattern in the spurs 
has been recently discussed by Vidal et al. (2015). Here we wish 
to draw attention to the regions of organized field orientation at 
high latitude, outside the NPS, which parallels it quite closely 
over most of its length. Some of this region is occupied by the 
cold border, where the magnetic field may be organized by the 
expansion of the loop, e.g., as modelled by Vidal et al. (2015). 
However, Fig. 24 shows that the parallel-field region extends 
past longitudes of 60°, well beyond the outer limit of the cold 
border. 

The origin of this high-latitude, inter-loop polarization is un¬ 
clear. Wolleben (2007) sees it as emission from a “new loop”, in 
whose shell the solar system is embedded, but the evidence for 
this loop is weak (Sect. 5.4) and its proposed interpenetration of 
Loop I is, to say the least, dynamically problematic. In a some¬ 
what similar picture, Vidal et al. (2015), following Heiles (1998), 
place the Solar System on the surface of Loop I, so that loop 
emission fills a hemisphere: this can explain the high-latitude 
polarization, but it is not clear why the cold border seems to 
have a significantly smaller radius than 90°. If the loop is nearby 
and the polarization is from the halo, the alignment with Loop I 
would be an accident. However, especially if the loop is as large 
as suggested in Sect. 5.3.2, the emission may be from the imme¬ 
diate environment of the loop, suggesting that Loop I is bright¬ 
est where the interstellar field is in the plane of the shell surface. 
This is exactly as predicted by Spoelstra (1972), who applied the 
model of van der Laan (1962) to the NPS. This appeals to the 
asymmetry of magnetic stress, which has a net tension along the 
field lines, so magnetic forces oppose expansion only perpen¬ 
dicular to the field. The same model has been invoked to explain 
the occurrence of barrel-shaped SNRs aligned with the Galactic 
plane (Gaensler 1998). 

Much of the random component of the Galactic magnetic 
field is on scales smaller than Loop I (e.g., Haverkorn et al. 2008; 
Brown & Taylor 2001), so the external field around it is likely 
to show significant variation in magnitude and direction. Such 
large objects are therefore unlikely to show the bilateral sym¬ 
metry of classic barrel-shaped remnants. Nevertheless the basic 
physics still operates; thus, segments of the shell expanding per¬ 
pendicular to the field, along the edge of the NPS, will be more 
strongly impeded and so require higher internal pressure, and 
therefore stronger nonthermal and thermal (X-ray) emissivity. 
By the same token, we expect the fainter parts of the shell to 
expand fastest, consistent with our interpretation of the bottom- 
right arcs as a bulge in the loop boundary. Of course, asymmetric 
densities in the ISM will also play a role, and as noted above, if 
the loop is at D 140 pc the observed density gradient is in the 
required sense. 

In old supernova remnants the projected magnetic field is 
commonly aligned with the outer surface, i.e., the shock propa¬ 
gating through the ISM. This alignment is naturally produced by 
shock compression, which ensures that the field lines behind the 
shock tend to be parallel to the shock plane. In this picture there 
is no global field order in the shock plane, but the field anisotropy 
causes polarization at viewing angles other than face-on, reach¬ 
ing a maximum when the line-of-sight is in the shock plane (see 


31 











Planck Collaboration: Planck diffuse low-frequency Galactic foregrounds 


Laing 1980). Applied to Loop I, this model implies that both 
the NPS and the internal spurs are tangential views of shock 
fronts. There is no problem with this if the internal spurs are 
background features, as suggested for one of the most prominent 
in Sect. 5.5 below. However, if, as seems likely, at least some of 
the internal spurs are physically associated with Loop I, then 
we have to explain multiple concentric shock waves. The ex¬ 
pected supernova rate for the Sco-Cen OB association is around 
1 Myr“^ (de Geus 1992), while remnant lifetimes are conven¬ 
tionally < 30kyr, and in the low-density environment of a su¬ 
perbubble they are likely to be shorter. Therefore multiple ac¬ 
tive supernova remants within the superbubble are very unlikely, 
even if Loop I is powered by a much more massive OB associ¬ 
ation than Sco-Cen. However, if a re-energizing blast wave has 
recently run into the contact discontinuity at the outer surface of 
the superbubble, it will separate into a transmitted and reflected 
shock, and for any non-ideal geometry a rather complex pattern 
of crossing shocks is expected in the reflected wave. MHD sim¬ 
ulations of this scenario would be illuminating, but are beyond 
the scope of this paper. 

A very different model for the NPS has been proposed by 
Heiles (1998) and is implicitly applied to the internal filaments 
of Loop I by Vidal et al. (2015). In their geometry the NPS is 
not a tangential view of the superbubble boundary, since the so¬ 
lar system is on the bubble surface. The organized held pattern 
is due to the wrapping of a relatively ordered interstellar held 
over the surface of the expanding bubble. The spurs are bundles 
of held lines with an enhanced density of cosmic ray electrons 
and therefore enhanced synchrotron emission. An attractive fea¬ 
ture of this is that a very simple geometric model gives a good 
qualitative fit to the held pattern in the loop. However, the model 
has the peculiarity that the emission is implied to be dominated 
by the far side of the loop; the near side, in which the viewer 
is embedded, would contribute a parallel polarization over the 
whole hemisphere, quite unlike the observed pattern. We note 
that in this model, Sh 2-27 can depolarize the loop emission even 
though located near the centre of the bubble. 

A hybrid of the two models is possible if the solar system 
is moved off the loop surface in the Heiles geometry, so that 
the NPS resumes its usual role as the projected loop boundary. 
However, this would worsen the agreement of the model and 
observed held pattern, and is still subject to the problem that 
emission from the near surface of the bubble would give a much 
more uniform projected held. 

5.3.4. Possible contamination of the CMB by Loop I? 

Liu et al. (2014) have recently argued that emission from 
Loop I is contaminating current microwave background maps. 
Specifically, they demonstrate an alignment between the BHS 
small circle and positive peaks in the low-multipole {i < 20) 
structure of the WMAP ILC map (Bennett et al. 2013a); von 
Hausegger et al. (2016) confirm that the alignment is also seen 
in the Planck CMB maps.^^ showed that the analysis by Liu et 
al. underestimated the the probability, p, of a chance alignment, 
but von Hausegger et al. (2016) present a stronger statistical test 
that gives p < 3 x 10“^. The amplitude of these peaks is about 
ten times larger than worst-case errors due to foreground con¬ 
tamination in the ILC maps estimated by Bennett et al. (2013a), 
which occur around the edge of the mask recommended for use 
with their ILC map at low multipoles (Liu et al. applied no mask 
in their analysis). The most convincing alignment with the cir- 

Ogbum(2014) 


cle is in the southern Galactic hemisphere, including along the 
bottom-right section that we note above shows no sign of emis¬ 
sion from the loop border. In general these peaks in Liu et al’s 
low-/ map show little correlation with the actual synchrotron or 
dust emission from Loop I, especially if we mask the Galactic 
plane as Bennett et al. recommend. We note that Liu et al. (2014) 
suggest that the contamination might be due to magnetic dipole 
emission from dust grains associated with Loop 1. However, the 
securely-detected dust associated with Loop I is mostly located 
at larger radii than the synchrotron ridge, and would not itself 
give a significant signal in the analysis of Liu et al. (2014) or 
von Hausegger et al. (2016). 

Moreover it seems unlikely that the dust in Loop I is par¬ 
ticularly unusual; even if we posit a unique dust type in the 
loops, Mertsch & Sarkar (2013) have argued convincingly that 
similar structures are scattered throughout the Galactic disc and 
contribute a significant fraction of its synchrotron emission. By 
the same token, the proposed dust emission from distant loops 
should accumulate along the Galactic plane, but no such signal 
is apparent in the CMB maps, von Hausegger et al. (2016) ar¬ 
gue that this expected signal is restricted to \b\ < 5°, a region 
where component separation is unreliable (hence the need for a 
Galactic mask). While it is true that it is not possible to accu¬ 
rately estimate the CMB signal in this region, we can rule out 
foreground contamination at the level required by this model. 
The integrated brightness of the distant loops is is an order of 
magnitude larger than that of Loop I, and the features responsi¬ 
ble for the Loop I alignment are among the brightest in the CMB 
maps after Altering to I < 20, with amplitudes 50-100yuK. Using 
this to scale the arbitrary units in Figure 5 of von Hausegger et 
al., the brightness from distant loops on the Galactic plane aver¬ 
aged over |/| < 60° should be 900yuK, and is still over 100/tK 
at \h\ = 4°. Fig. 26 shows profiles of derived CMB brightness 
against b for the four main Planck component separation meth¬ 
ods. Except for SEVEN, which is the method least well adapted 
to removing diffuse foregrounds, substantial contamination from 
Galactic foregrounds is confined to \b\ < 2°. Even the SEVEN 
map has smaller Galactic residuals than predicted from the von 
Hausegger et al. analysis. 

In short, there is little reason to believe in the contamination 
suggested by Liu et al. The alignment they suggest is with a theo¬ 
retical schematic rather than the observed synchrotron loop, and 
is even further away from the actual dust emission in the loop. 
The formal improbability that the observed alignment occurs by 
chance, p ^ 0.2% in the best (i.e., Planck) CMB maps, should be 
offset against the low prior probability that this schematic acci¬ 
dentally aligns much better with a previously-undetected struc¬ 
ture than with the synchrotron loop that it was intended to model. 
This p-value does not allow for any look-elsewhere effect such 
as potential alignments with other loops. Moreover, their pro¬ 
posed physical model implies additional strong contamination 
along the Galactic plane, which is not observed. 

5.4. Other loops and spurs 

Three other synchrotron loops that are generally agreed upon 
(e.g., Vidal et al. 2015), and there are a dozen or so shorter fil¬ 
aments or spurs, with several mutually-incompatible proposals 
for joining some of them into loops. Given our emphasis on the 
non-circularity of Loop I, and the fact that it seems all too easy to 
And small circles passing plausibly close to a sequence of ridges 
on the sky maps, we do not And any of these newly proposed 
loops convincing; however, their component spurs are certainly 
real. We discuss the more interesting examples in turn. 
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Galactic latitude [°] 

Fig. 26. Mean brightness vs. Galactic latitude, b, in the Planck 
CMB maps from Planck Collaboration IX (2016), in the range 
|/| < 60°. Solid green: NILC; dotted blue: SMICA; dot-dashed 
magenta: SEVEN (all at 5 arcmin resolution); dashed red: base¬ 
line Commander model at 1° resolution (Planck Collaboration X 
2016). The first three methods all use inpainting to fill masked 
regions around bright compact features in some or all of the in¬ 
put maps: pixels affected by this have been omitted from the 
averaging. High amplitudes close to /? = 0 are due to residual 
foreground contamination. 


Loop II. Also known as the Cetus arc. Loop II is a very dif¬ 
fuse, 45°-radius structure centred at (/, Z?) = (100°, -32?5). The 
faint, broad emission on the left-hand side of the maps in Fig. 2 
is attributed to it. As originally defined by Large et al. (1962), 
the right-hand end of the arc was considered to be the “South 
Polar Spur” (SPS), which descends from the plane at / = 45° 
(visible in Figs. 20 and 21). The SPS is prominent in both the 
unsharp masked 408 MHz map and the synchrotron polarization 
map, but curves towards the right as it leaves the plane (as does 
its projected magnetic field), in the opposite direction from the 
notional path of Loop II; it therefore appears to be unrelated. The 
BHS small circle for Loop II just crosses the plane, but no emis¬ 
sion from it has been traced in the northern Galactic hemisphere; 
the top of the loop may be obscured by Galactic plane emission, 
but it is also possible that the expanding shell was halted by the 
dense gas in the plane. Alternatively, Weaver (1979) and Heiles 
(1998) suggest that Loop III is its northern extension, although 
the two loops are offset by 24° in longitude. No reliable distance 
information is available for Loop 11. 

Figure 24 shows a coherent field pattern roughly parallel to 
the locus of Loop II that is not visible in the WMAP data. The 
polarized intensity essentially disappears into the noise at the 
bottom of the loop near the South Galactic Pole (SGP), but a 
coherent polarization field can be traced along a path that goes a 
few degrees further towards the SGP than the BHS fit, and closes 
a few degrees higher in longitude from the SPS. Some of this 
polarization is detected at 1.4 GHz by Wolleben et al. (2006), 
but the polarization angles are more disordered, presumably due 
to Faraday rotation, and so the alignment with the loop is not so 
clear. 

Loop III. Loop III is centred at (/,/?) ^ (124°, 15?5); BHS give 
a radius of 32?5, although Figs. 20 and 24 show that polar¬ 
ized emission is detected several degrees further out. Loop III 



Fig. 27. Combined polarization map of the Cygnus spur, the 
right-hand end of Loop III. The map is centred at (/, Z?) = 
(90°, +21°) and the graticule has a 15° spacing. 


is clearly detected around most of its perimeter in the north¬ 
ern hemisphere in both total intensity and polarization. Its right- 
hand spur, rising from the plane directly above the Cyg X re¬ 
gion, passes through the region of deep coverage near the North 
Ecliptic Pole (/, b) ^ (96°, 30°) on the Planck maps, giving a par¬ 
ticularly clean view of the well ordered magnetic field (Fig. 27). 
This Cygnus spur is relatively narrow (FWHM 5°) until it 
reaches ^ 30° after which the loop becomes much broader and 

fainter. As it returns to the plane its polarization is obscured by 
the strongly polarized Fan region (100° < I < 170°), in which 
the magnetic field is parallel to the plane and therefore orthog¬ 
onal to the expected field in the loop. BHS argue that the loop 
re-emerges from the fan region in the southern hemisphere at 
I ^ 150°, and the 5-field pattern in our synchrotron polarization 
map is consistent with this, but the loop does not convincingly 
close in the south. Vidal et al. (2015) imply that their filament 
Ills may be associated with Loop III proper, but this is one of 
the regions where WMAP and Planck disagree the most, and the 
existence of this filament is not confirmed in the Planck maps, 
except for a short section near (/, b) = (82°, -30°) (see Fig. 24). 
Kun (2007) discusses possible interactions between Loop III and 
the cold gas north of the plane, notably the North Celestial Loop 
(Meyerdierks et al. 1991) that may be swept-up material. 

Loop IV. Loop IV is centred at (/, b) = (315°, 48.5°) and is pro¬ 
jected entirely inside Loop 1. Only an arc subtending 20° along 
its low-latitude rim is clear in our polarization map; see Vidal 
et al. (2015) for analysis. The high-latitude rim of Loop IV that 
parallels the top of the NPS is very noisy here but clear in the 
21cm polarization map of Wolleben et al. (2006). 

The South Polar Spur (SPS) and Filament X. As discussed 
above, the SPS was once seen as the low-longitude edge of Loop 
II, but in fact curves in the opposite direction so the outside edge 
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Fig. 28. The South Polar Spur. Left: P/a^c^-WMAP polariza¬ 
tion map; right: Planck 353-GHz dust polarization. The maps 
are centred at (/,/?) = (50°,-45°) and the graticule has a 10° 
spacing. 


of the SPS is to its right. It also forms the most prominent seg¬ 
ment of the SI loop proposed by Wolleben (2007), and was la¬ 
belled filament Vllb by Vidal et al. (2015). Wolleben’s view of 
the SPS was obscured by Faraday rotation near the plane; as 
noted by Vidal et al. (2015), the high frequency (> 20 GHz) po¬ 
larization structure of the SPS is more strongly curved near the 
plane than expected for Wolleben’s path. The other section iden¬ 
tified by Wolleben as part of his new loop is usually seen as one 
of the interior ridges of Loop I and the northern part of Loop 
IV; thus, all components of loop SI appear to be parts of smaller 
structures and there is no clear evidence that SI exists as a co¬ 
herent physical structure. 

The SPS has its own cold border just outside the radio ridge, 
at I ^ 60°, -30° > b > -70°, which is particularly clear in dust 
polarization (Fig. 28). The border is visible in Hi at 0 < vlsr < 
10kms“^ (Fig. 21d,e), which implies a distance of less than a 
few hundred parsecs. 

Given our discussion above, we are reluctant to define a 
loop by extending a small circle from the clearly-detected arc 
of the SPS. In fact, a slightly fainter spur immediately north of 
the plane (Vidal et al.’s filament X) seems to be a reflection- 
symmetric twin of the SPS, suggesting a wasp-waisted cavity 
shaped by the dense gas in the plane (see the left-hand edge of 
Fig. 21c). 

Our map allows filament X to be traced further north, to 
b = 37°, where it curves over and enters the region of bright 
diffuse polarized emission just outside the NPS (Fig. 25). 
Unusually, the field in this northern section is perpendicular to 
the filament. This alignment is not due to a contribution from 
the underlying diffuse emission; if it were, the filament would 
appear as a trough (due to cancellation), not a ridge in polarized 
intensity. A possible return section is marked in Fig. 20. 


The Orion-Taurus Ridge and Filament XI. Unlike the other 
spurs, these features do not emerge from the Galactic plane but 
run parallel to it at closest approach, at (/,/?) (220°,-15°), 

and (/, b) ^ (230°, 40°), respectively (see Fig. 20). They are vis¬ 
ible in both unsharp-masked total intensity and in polarization 
(see figure 2 of Vidal et al. 2015). Milogradov-Turin & Urosevic 
(1997) and Borka et al. (2008) propose the existence of a pair of 
very large, overlapping loops (V and VI in their designation) that 
connect a sparse collection of features including the right-hand 
end (at I < 180°) of the Ori-Tau ridge and Vidal et al.’s filament 
Ills. However, the rest of the Ori-Tau ridge curves far away from 
the locus of Loops V/VI, and the magnetic field continuously 
follows that curvature, so these loops appear to be spurious. 

Filament XI is an approximately 80°-radius arc that seems 
to curve away from the plane at both ends (Fig. 20); however, 
it is unclear if features near I = 264° belong to the left-hand 
end of Filament XI or the right-hand side of Loop I (or filament 
XII, if not part of Loop I). Extended to a full loop. Filament XI 
encloses most of Galactic quadrant 3 north of the Plane, which 
is one of the faintest regions of the polarized sky; over much of 
this quadrant no polarization is clearly detected by Planck and/or 
WMAP, even at 4° resolution. 

/ = 45° feature. This is the highly-polarized patch covering 
12° < /? < 45° and from the edge of the NPS to I ^ 50° (Fig. 25). 
The low-/? edge of the polarization feature coincides with a rapid 
brightening of the inner halo in total intensity, and when polar¬ 
ization re-appears at /? ^^ 8° it has an orthogonal orientation, 
roughly parallel to the Galactic plane. Hence this minimum is 
due to a cancellation between the plane-parallel field in the in¬ 
ner halo and the near-vertical field in our feature, and does not 
mark the physical edge. A faint ridge runs through the top half 
of the feature at / 45°, which is just visible in the 408 MHz 

map (Fig. 21a,b); otherwise there is little trace in total intensity. 
The magnetic field runs at a slight angle to the ridge of the NPS 
throughout; in the adjacent section of the NPS the projected field 
is misaligned with the ridge-line and in better agreement with 
that in the / = 45° feature. Plausibly the feature overlaps in pro¬ 
jection with Loop I, but if so the current data do not allow us to 
trace it inside the loop. 

Smaller loops. Mertsch & Sarkar (2013) have emphasized that 
analogues to the large synchrotron loops should be common, and 
more distant and smaller examples should therefore be visible in 
the sky maps close to the Galactic plane. One rather clear ex¬ 
ample is the polarization structure south of the CygX region, 
centred at (/, b) = (90°, -5°), which seems to be a superposition 
of a 5°-radius loop around this bright star-forming complex and 
another spur heading to the lower right (Fig. 29). The radius is 
about 120pc, assuming a distance of 1.4kpc for CygX (Rygl 
et al. 2012). A northern counterpart is less obvious, partly be¬ 
cause of the two bright supernova remnants, HB 21 and W 63, 
that lie on the likely path. Figure 29 also shows another possible 
distant loop: the set of three arches in the southern hemisphere, 
centred near I = 60°; they are also marked in Fig. 20. The coher¬ 
ence of the polarization vectors allows the outer arch to be fol¬ 
lowed around to (/, b) = (17?4, -12?9), just south of the Cygnus 
Loop supernova remnant. While the polarized intensity suggests 
that the inner arch returns to the plane at / 61°, the polariza¬ 

tion vectors are orthogonal to the apparent ridge-line, and so this 
might be a different structure. Figure 29 also shows that the dust 
polarization has the same magnetic pattern as the synchrotron 
emission around Cyg X. There is also a polarized dust filament 
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Fig. 29. Top: combined polarization map with magnetic field 
vectors overlaid, centred on (/, Z?) = (70°,-7?5) showing sev¬ 
eral polarized arches south of the plane. Bright compact fea¬ 
tures along the plane are all SNRs: DA 530 (G93.3-f6.9); HB 21 
(G89.2-r4.7); W63 (G82.2-r5.3); CTB 87 (G74.9-fl.2); the 
Cygnus Loop (G74.0-8.5); CTB 80 (G69.0-f2.7); and W51C 
(G49.1-0.6). The large spur at I ^ 50° is the SPS. Bottom: dust 
polarization in the same region as seen at 353 GHz. Maps are 
rectangular projections and have an asinh colour scale. 


that runs roughly along the inner edge of the middle I = 60° 
arch, while the outer arch parallels the dust in polarization as it 
runs behind the SPS. 


while searching for a y-ray counterpart to the microwave haze 
(e.g., Finkbeiner 2004; Planck Collaboration Int. IX 2013) in the 
Fermi data. 

The origin of the bubbles is not clear, although their location 
and symmetry with respect to the Galactic centre suggest that 
they originate there. They have a y-ray spectrum significantly 
harder than the inverse Comp tom emission from the Galactic 
Halo or the one from pion decays from collisions of CR, with 
the ISM protons and heavier nuclei (Su et al. 2010). Different 
models have been proposed to explain their origin, most of them 
relating to a recent AGN-type activity at the Galactic centre. For 
a recent overview on the bubbles and their spectral behaviour in 
y-ray, see Dobler (2012) and Ackermann et al. (2014). 

We compare our combined Planck-y^MAP polarization am¬ 
plitude map at 28.4 GHz with a full-sky Fermi map from 
Ackermann et al. (2014). The y-ray map was produced using 
50 months of Fermi LAT data (Atwood et al. 2009), using only 
the “UltraClean” class events, i.e., the sample of events with less 
contamination from misclassified interactions in the Fermi LAT 
instrument. See Ackermann et al. (2014) for a detailed descrip¬ 
tion on the Fermi data analysis. We use the high energy map, 
which covers 10-500 GeV. In this energy range, the emission 
from the Fermi bubbles appears clearly without any component 
separation at latitudes h > 20°. The main foreground emission 
in gamma rays is from decay of particles produced by CR 
protons interacting with the ISM (Ackermann et al. 2012). If the 
CR density and spectrum are roughly spatially uniform over the 
Galaxy, then the emission will be proportional to the ISM 
column density (Su et al. 2010). A good tracer for the ISM col¬ 
umn density are maps of thermal emission from dust grains, 
since dust is well mixed in the ISM and its emission is opti¬ 
cally thin. Here, we use the Planck 353 GHz optical depth map 
from Planck Collaboration XI (2014) as a column-density proxy 
for 71^ emission to fit and remove it from the y-ray map. 

In the middle panel of Fig. 30, we show the resulting 10- 
500 GeV Fermi map. A region of about +5°along the Galactic 
plane, where the subtraction of the scaled Planck 353 GHz map 
left large residuals, has been masked out and is shown in grey. 
The left panel in Fig. 30 shows our combined P/a^c^-WMAP 
polarization amplitude map, where narrow filaments are visi¬ 
ble that correspond to the border of the Fermi bubbles. We note 
that the northern filament is much clearer than the southern one. 
The ridge-lines of these filaments is over-plotted in the Fermi 
map on the right to show the correspondence. It is remarkable to 
see how closely the polarized filament follows the border of the 
Fermi bubbles. We also show on the right panel the denoised 
and source-subtracted Fermi bubbles template constructed by 
Selig et al. (2015) using a Bayesian inference algorithm. Their 
reconstruction has a poorer angular resolution than the less pro¬ 
cessed version we show at the centre of Fig. 30, but still it is 
clear that the synchrotron filaments outline the border of the 
bubbles. The polarized filaments are unresolved in the Planck 
polarization maps (40' FWHM beam), implying that these are 
very narrow synchrotron structures. This is further confirmation 
of a sharp border for the Fermi bubbles, and models that fail to 
predict this feature are now less favoured (see Dobler 2012 for a 
description of some of these models). 


5.5. Fermi bubbles in polarization 

The Fermi bubbles are two large structures extending perpen¬ 
dicular to the Galactic plane, up to about +55°. They were 
first discovered by Dobler et al. (2010) using spatial templates 


Polarization spectral index. We have measured the polarization 
spectral index, y5poi, between the bias-corrected WMAP 23 GHz 
and Planck 30 GHz using T-T plots of a region that encom¬ 
passes the north filament. The left panel of Fig. 30 shows two 
regions in the polarization intensity map that we selected to mea- 
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Fig. 30. Left: cartographic projection of the combined Planck-^MP^ polarization amplitude, centred at the Galactic centre, with 
graticule lines at 15° intervals. The regions defined by red and yellow dashed lines were selected to calculate a polarized spectral 
index; in red is the filament around the Fermi bubble and in yellow a control area. Centre: cartographic projection of the 10- 
500 GeV map from Ackermann et al. (2014), with the emission subtracted as described in the text, showing the Fermi bubbles. 
Right: Fermi bubbles component from Selig et al. (2015). The black outline corresponds to the centre of the narrow filaments visible 
in the polarization map on the left. 


sure ySpoi, one including most of the northern filament (in red) and 
a second region that we use as control (yellow). 

The T-T plots are shown in Fig. 31. The polarization spec¬ 
tral index of the Fermi bubble filament is ySpoi = -2.54 + 0.16, 
while the nearby control region has ySpoi = -2.90 ± 0.36. The 
spectral index of the region that includes the filament is much 
flatter than the standard value of -3.0 normally assumed for dif¬ 
fuse synchrotron emission. It is also flatter than all the values 
presented in Fuskeland et al. (2014) and Vidal et al. (2015), who 
explored the polarized spectral indices in different regions us¬ 
ing WMAP data. This flatter spectral index indicates that the 
energy distribution of the synchrotron radiation along the fil¬ 
ament has more higher energy electrons than that of the dif¬ 
fuse component next to it. The polarization spectral index of 
the Fermi bubble filament we measure here is identical to the 
spectral index of the microwave haze of y^Haze = 2.54 ± 0.05 
(Planck Collaboration Int. IX 2013). This supports the relation¬ 
ship between the filament and the haze/bubble emission. 


Fermi bubbles and Loop I. A relationship between Loop I (see 
Sect. 5.3.1) and the Fermi bubbles has been hypothesized (e.g., 
Kataoka et al. 2013). This would put Loop I at the Galactic cen¬ 
tre, with a much larger size and radio luminosity, as long argued 
by Sofue (e.g., Sofue 1977, 1994, 2015). 


With the new Planck-WMAF combined polarization map, 
in Sect. 5.3.1 we trace Loop I below the Galactic plane. It is 
clear from the top panel of Fig. 20 and from Fig. 30 that the 
southern part of the Fermi bubble extends outside Loop I, and 
there is no trace of an interaction with the bubble in the radio 
maps, given the continuity and smoothness of the Loop I polar¬ 
ized filaments in the southern hemisphere. The two structures 
must therefore be at different distances. The Fermi bubbles are 
well centred on the Galactic centre, show a pinched structure 
symmetric about the Galactic plane, as expected for an outburst 
from the nucleus (Sofue 1994), and are unique in the y-ray sky. 
All this makes it highly likely that they are located at the distance 
of the Galactic centre, as usually assumed. In contrast. Loop I is 
centred 35° from the nucleus and significantly above the plane, 
while on its left-hand side, where it can be clearly followed, it 
appears to extend through the plane without any sign of devia¬ 
tion (Sect. 5.3.1). This is just as expected if it is relatively nearby 
and embedded in the disc rather than extending far beyond it. Of 
course, it also has a larger angular size than the Fermi bubbles, 
and there are at least several similar loops (Sect. 5.4). We are 
therefore confident that Loop I is a foreground feature. 
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Filament Fermi Bubble 




WMAP 22.8 GHz P [/zK] 

Fig. 31. T-T plots between the polarization amplitude maps of 
WMAP K-band and Planck 28.4 GHz. The two plots correspond 
to the regions defined in the left panel of Fig. 30. The Fermi 
bubble filament shows a flatter spectrum than the control region. 




-2.3 -1.0 log (mK) -2.8 ^ -1.3 log (mK) 


Fig. 32. Top left: Ha map in the velocity range -80 < Visr < 
-40kms“F Notice the vertical filament that runs at I ^ 75° 
and -60° < b < -20°. Top right: Faraday depth map 

from Oppermann et al. (2012), which also shows a filament 
at the same location, with a mean value of -25 rad m~^ along 
its extension. The filament has a counterpart in “absorption,” 
visible as a trough in polarization intensity maps. Bottom 
left: WMAP 23 GHz polarization intensity map. Bottom right: 
Planck-WMAF polarization intensity map. The Her map on the 
Top left comer has an angular resolution of 1°, while the other 
three maps have a common resolution of 3°. The grid spacing is 
15? 


5.6. An anti-correlation of an Ha filament and polarized 
intensity 

Another interesting aspect that can help in the determination of 
the distance to some of the polarized structures is the discov¬ 
ery of a region that shows an anti-correlation between Her emis¬ 
sion and the polarized intensity in our combined Planck-WMAP 
map. This has not been noticed before as far as we are aware, 
because the filaments are located at high Galactic latitude, in the 
Faraday-thin regime. 

In Fig. 32 we show an Her map, integrated over the velocity 
range -80 < Visr < -40kms“\ from the Wisconsin H-Alpha 
Mapper (WHAM) survey (Haffner et al. 2003). The filament that 
runs along I ^ 75° is about 40° in length. The WMAP 22.8 GHz 
polarized intensity, and the Planck-^MAP polarized intensity 
map also shown on Fig. 32, show a trough at the same location as 
the Her filament. This feature is also visible in the Faraday depth 
map from Oppermann et al. (2012), shown on the same figure. 
The Faraday depth, f is proportional to the magnetic field along 


the line of sight and the electron density: 

(p = kF ^ drn^{r)Br{r). (5) 

Jro 

We have ruled out a chance correlation since the polarized 
feature is also clearly visible, at higher angular resolution, in the 
maps from the GALFACTS (Taylor & Salter 2010) consortium 
(Jereon Stil, priv. comm.), and the correlation is still strong at 
angular scales of a few aremin. 

The origin of the observed anti-correlation between Ha in¬ 
tensity and polarization amplitude is not clear. A first possibility 
is that the trough visible in the combined Planck-WMAF po¬ 
larization map is a depolarized region, meaning that the Ha fil¬ 
ament lies in between the polarized background emission and 
us. If this is the case, the ionized gas traced by the Ha map 
could produce Faraday rotation, due to the presence of a mag¬ 
netic field in the plasma, depolarizing the diffuse background 
emission along its extension. However, this hypothesis is not 
compatible with the low density of the ionized gas. The inten¬ 
sity of the filament in Ha is 2 R above the background at the 
original resolution (6') of the Finkbeiner (2003) Ha map. This 
corresponds to a mean electron density of 2.0 cm“^ (assuming 
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that the filament has a thickness of 1 pc), which is very low to 
produce significant Faraday rotation at 23 GHz for typical val¬ 
ues of the magnetic field. Moreover, the Faraday depth map from 
Oppermann et al. (2012) has a value of about 25radm“^ along 
the filament, which corresponds to a 0?3 change in polarization 
angle at 23 GHz. Therefore, Faraday rotation of this high lati¬ 
tude filament is not enough to cause any major depolarization at 
23 GHz. 

A second alternative is that there could be a strong coherent 
magnetic field parallel to the line of sight along the filament. 
This might be the result of the dynamical processes that created 
the filament. If this is the case, there would be less polarized 
emission at the filament location than around it. 

Lastly, it might well be that the synchrotron emission from 
the region of the Ha filament is intrinsically weakly polarized. 
This might be due to a less organized magnetic field in this re¬ 
gion in comparison with the diffuse emission seen in the vicinity 
of the filament on the polarization map. 

We also note the fact that the Ha filament is only visible at 
negative radial velocities. This corresponds to the velocity range 
of the Perseus arm of the Galaxy. If the Ha filament belongs 
to that arm, it would imply that the distance to the diffuse syn¬ 
chrotron background is much larger than a few hundred parsecs, 
lying at least 2kpc away from us. 

5.7. Limits on AME polarization 

If AME is solely due to spinning dust particles, then we ex¬ 
pect it to have a very low polarization percentage. The level of 
spinning dust polarization depends on the alignment efficiency 
of the small grains and PAHs in the interstellar magnetic field. 
Lazarian & Draine (2000) considered resonance paramagnetic 
relaxation, which predicts < 1.5 % polarization for frequencies 
> 20 GHz; Hoang et al. (2013) use constraints on the alignment 
of grains seen in ultraviolet polarization to predict that AME 
will have a polarization of < 0.9 % at frequencies above 20 GHz. 
At lower frequencies (< 10 GHz), the polarization fraction will 
be higher, but as the AME spectrum steeply decreases at low 
frequencies while other foreground components are increasing, 
this will be more difficult to detect than at the peak of the AME 
spectrum. Given this, and the fact that the observed polarization 
percentage would naturally be less than this due to beam and 
line-of-sight depolarization, we do not expect to detect signif¬ 
icant AME polarization with Planck, i.e., the polarization per¬ 
centage should be < 1 %. 

Observational constraints on AME polarization have so far 
been placed using relatively compact, isolated clouds, where 
AME is known to be strong. Rubino-Martm et al. (2012a) review 
the available constraints and their relation to theoretical mod¬ 
els. The best constraint comes from the Perseus region, which 
has significant AME emission and relatively little contaminat¬ 
ing synchrotron emission. Battistelli et al. (2006) reported a 
weak detection in Perseus at 3.4^}| % at 11 GHz, while later 
measurements obtained by Lopez-Caraballo et al. (2011) using 
WMAP data found a 2cr limit of < 1 %; Dickinson et al. (201 1) 
have also measured a Icr limit of <1.4% for Perseus (as well 
as an upper limit of <1.7% in p Ophiuchus), and Genova- 
Santos et al. (2015) found a Icr limit of <2.8 % in Perseus at 
19 GHz. Recently, Battistelli et al. (2015) have claimed to de¬ 
tect polarized AME emission at 21.5 GHz from RCW 175, an 
Hii region, where they measured a polarization percentage of 
(2.2 + 0.2 (random) + 0.3 (systematic)) % at 21.5 GHz. However, 
Battistelli et al. (2015) argue that a large fraction of this sig¬ 
nal could be residual synchrotron radiation, leaving little or no 


Perseus 


Plume 




Fig. 33. Top: Coimnander synchrotron component maps; mid¬ 
dle: Commander AME component maps; and bottom: debiased 
weighted polarization maps. Left: Perseus; right: the Pegasus 
plume. The colour scales are linear here. The regions are defined 
and masked as per Sect. 4.2. The Perseus region does not contain 
significant polarized synchrotron emission, while the Pegasus 
plume contains several polarized synchrotron arcs. 


polarized AME. There have been very few attempts to con¬ 
strain the polarization from the large-scale diffuse AME; Kogut 
et al. (2007) find that AME accounts for less than 1 % of the 
total polarized signal, while Macellari et al. (2011) placed a 
limit of < 5 % based on template fitting of WMAP data. Planck 
Collaboration Int. XXII (2015) placed upper limits by corre¬ 
lating WMAP and Planck polarization maps with the Planck 
353 GHz dust polarization map. They found that the correlation 
turned upwards at low frequencies, and at WMAP K-band they 
placed an upper limit on AME polarization of 16%; they also 
noted that the upturn could be explained by dust-correlated syn¬ 
chrotron polarization rather than AME polarization. 

We use the debiased polarization map assembled using a 
weighted average of both Planck and WMAP data (see Sect. 5.1) 
to look for polarized emission that is correlated with diffuse 
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Fig. 34. T-T plots between the polarization and AME intensities 
(both in units of p,K at 22.8 GHz). Top: Perseus; bottom: Pegasus 
plume. The scatter is mostly due to thermal noise in the polar¬ 
ization map (uncertainties in the polarization intensity in these 
regions are typically around 5 p,K). The fitted slope gives the po¬ 
larization fraction for the region, based on the linear fit (dashed 
line). 


AME from the Commander solution. We note that the weighted 
polarized map assumes a spectral index of (3 = -3.0, rather than 
an AME-like spectrum. We create T-T plots between the AME 
and polarization intensities, both rescaled to p,K at 22.8 GHz, 
in order to measure the percentage polarization. Assuming that 
any potential AME polarization has a constant percentage across 
the region, then this should result in a linear correlation of the 
data points in the T-T plots. We use the uncertainty maps and 
assume a model uncertainty of 10 % for both the Commander 
AME solution and the weighted polarization map; the latter un¬ 
certainty is subdominant in these regions. We fit for the slope of 
the correlation, which directly gives the percentage polarization. 
We repixelize the maps to Aside = 64 to ensure that the pixels are 
independent. We note that the debiasing method used will not 
perfectly remove all of the noise bias, particularly in regions of 
low S/N, but any residual bias should be lower than the thermal 
noise limits where the S/N ratio is low. 


We focus on two high latitude regions from Sect. 4.2, namely 
Perseus and the Pegasus plume. We show the AME and debi- 
ased polarization maps in Eig. 33, and the AME-polarization 
T-T plot for each region in Eig. 34. AME dominates the to¬ 
tal emission at around 20-30 GHz in both of these regions. 
Unlike regions such as Corona Australis and Musca, which are 
strongly contaminated by highly polarized Galactic synchrotron 
structures (see Sect. 5.2), Perseus and the Pegasus plume are 
in parts of the sky with relatively little synchrotron polariza¬ 
tion, although the Pegasus plume does have two polarized syn¬ 
chrotron features crossing the region; these can be seen in Pig. 33 
and are marked on Pig. 20. One of these synchrotron arcs is 
somewhat aligned with the AME plume, however it is unre¬ 
lated to the plume, since it extends beyond the plume towards 
the Galactic plane. The morphology of the polarized and total 
intensity synchrotron emission is quite different in the region of 
the plume, and the aligned arc is either highly polarized or has a 
fiat spectrum, since it is faint in the total intensity map. Unlike 
the Perseus and Pegasus plume regions, the Orion and T Orionis 
regions contain strong free-free emission: while free-free emis¬ 
sion is unpolarized there may be residual leakage from total in¬ 
tensity to polarization. In addition there is a broad background 
of polarized synchrotron emission in these regions, and due to 
their proximity to the Galactic plane there will be a correlation 
between synchrotron and thermal dust emission, simply because 
both components will be brighter towards the Galactic plane than 
they are off-plane. We apply the point source mask described in 
Sect. 4.2: in the regions of interest this only consists of masking 
one point source near the Pegasus plume. 

We find that the strongest constraint comes from Perseus, 
where we measure an AME polarization percentage (compared 
with the overall polarization intensity) of 0.6 + 0.5 %, which 
gives a 2(T limit of < 1.6 %, matching previous constraints. The 
Pearson correlation coefficient is 0.20. An improved Commander 
estimate of the AME in this region (see Sect. 4.2) would de¬ 
crease this upper limit by around 40%. In the Pegasus plume 
region we find a polarization percentage of 9.0 + 1.9 %, with a 
Pearson correlation coefficient of 0.29, which, given the pres¬ 
ence of correlated polarized synchrotron emission, we interpret 
to be a 2(T limit of < 12.8 %. We attempt to subtract some of 
the polarized synchrotron emission by assuming that the inten¬ 
sity is 10 % polarized, which removes some of the synchrotron 
emission from the region near the point source. This reduces the 
polarization percentage relative to the AME to 5.2 + 1.6%, or 
a 2 (t limit of < 8.4 %, with a Pearson correlation coefficient of 
0.24; assuming a higher level of synchrotron polarization results 
in significant negative pixels in the polarization map. We also 
compare the weighted low-frequency polarization map with the 
polarized dust emission at 353 GHz by separately comparing the 
Q and U maps at the two frequencies, at Aside = 64, in the re¬ 
gions of Perseus and the Pegasus plume. Similar analysis was 
performed by Planck Collaboration Int. XXII (2015) across the 
whole sky. We find similar results, with no significant detection 
in Perseus, and a weak detection in the Pegasus plume, which 
we ascribe to dust-correlated polarized synchrotron emission, as 
described above. 

In summary, we do not find strong evidence for diffuse polar¬ 
ized AME in Planck data, with a 2cr upper limit of 1.6 % 

in the Perseus region. In order to improve on current constraints, 
we will need either reduced noise levels in regions that are free of 
polarized synchrotron emission (such as Perseus) or the ability 
to subtract accurately the synchrotron emission from less clear 
regions. The forthcoming map of polarized synchrotron emis- 
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sion at 5 GHz by the C-Band All Sky Survey (King et al. 2014) 
will assist with the latter issue. 


6. Conclusions 

In this paper we have discussed the Galactic foreground emis¬ 
sion observed by Planck between 20 and 100 GHz. 

The total intensity emission comprises at least three distinct 
components: synchrotron; free-free; and AME, commonly at¬ 
tributed to spinning dust. These components emit smooth con- 
tinua that are typically decreasing with frequency. It is therefore 
extremely difficult to separate these components with any pre¬ 
cision. Furthermore, these individual components are often spa¬ 
tially correlated (e.g., towards the Galactic plane). 

We have applied an internal linear combination technique to 
Planck data to produce a map that is free of CMB, free- 
free, and thermal dust. This has given us some insight into the 
spectral form of emission at 23 GHz. For most of the intensity 
analysis we rely on the component-separation products from the 
Coinmander code. This involved fitting a parametric model to 
the data. However, a number of simplifications have had to be 
made, given the limited number of frequencies and complex¬ 
ity of the foreground emission. In particular, we have not so 
far been able to produce a detailed synchrotron map in inten¬ 
sity. Instead we have had to rely on a GALPROP spectral model 
to take into account the bulk of the synchrotron emission above 
20 GHz. In polarization the situation is somewhat clearer. The 
polarization maps at 20-50 GHz are dominated by synchrotron 
radiation; free-free and AME are expected to have very low po¬ 
larization fractions. 

Nevertheless, in this study we have been able to produce rea¬ 
sonable products for free-free and AME, from which we have 
inferred some physical results. Our polarization results do not 
rely directly on component separation. 

The main conclusions from our study are as follows. 

“ After applying an ILC procedure to remove CMB, free- 
free, and thermal dust emission at 22.8 GHz, we find that 
the residual emission has significant spectral variations. This 
emission is expected to be dominated by (weak) synchrotron 
and (stronger) AME emission. In particular, we find that for 
some regions the spectrum is flatter than that of free-free 
emission {{3 = - 2 . 1 ). 

- The aforementioned areas are largely towards bright Hii 
regions such as the f Ophiuchi cloud (Sh2-27)and the 
California nebula. Because synchrotron emission is not ob¬ 
served to be flatter than aboutyS = -2.3, the spectral changes 
are thought to be due to a high frequency component of 
spinning dust. Indeed, Planck Collaboration Int. XV (2014) 
found evidence for a spinning dust peak at around 50 GHz 
for the California nebula. 

“ Our free-free map appears to correspond well to previous 
maps made with WMAP data, and also to Her maps at high 
latitudes (where dust absorption is small). The amplitude rel¬ 
ative to Her is lower than the theoretical value for typical 
electron temperatures, as has been found in previous studies. 

- A natural explanation for the low free-free-to-Her ratio is that 
there is significant scattered Her light. We have attempted 
to estimate this fraction by correlating Her residuals (after 
subtracting the mean correlation) with the dust optical depth. 
We do indeed find a correlation with the thermal dust optical 
depth. If we assume that electron temperature varies by no 
more than +1000K, we find that the scattered Her emission 


accounts for 28 + 12%, or 36+12% after making a nominal 
correction for absorption by dust. 

- We also compared the Commander free-free map to the RRL 
Galactic plane survey of Alves et al. (2014). Again, we find 
a good overall correspondence between the maps. However, 
there are significant differences, which might be due to vari¬ 
ations in electron temperature. 

- The Commander AME map is closely correlated with ther¬ 
mal dust emission, as traced by FIR/sub-mm maps. We have 
identified several new diffuse regions away from the Galactic 
plane that appear to emit significant AME. These regions 
are typically associated with large dust cloud complexes. We 
have compared the emissivity of AME in these regions rela¬ 
tive to the amount of thermal dust emission. The AME emis¬ 
sivity against 545 GHz and T 353 varies by a factor of approx¬ 
imately 2 from region to region. The A Orionis region has 
a particularly high emissivity that warrants further investi¬ 
gation, ideally in conjunction with additional low-frequency 
data to reduce leakage between the free-free and AME com¬ 
ponents. 

- Comparison of the Commander foreground products with 
similar products from the WMAP team shows that the am¬ 
plitude of synchrotron emission is significantly lower in our 
analysis. This is particularly the case at low latitudes, where 
our fixed spectral model cannot account for the flatter in¬ 
dex reported by WMAP and in earlier Planck analyses. To 
compensate, our AME amplitude is significantly higher by 
a factor of around 2-4 (depending on which model is be¬ 
ing compared). The free-free component is in better agree¬ 
ment in general, but our model is typically lower by about 
10-30%, which has been noted before (Alves et al. 2010, 
2012). We emphasize that our model fits the data extremely 
accurately, with median residuals along the Galactic plane 
of < 1 % for Planck and < 2 % for WMAP channels (Planck 
Collaboration X 2016). 

- We have combined Planck polarization maps to con¬ 

struct a higher S/N ratio polarized intensity map. This is 
the highest S/N ratio synchrotron polarization map available 
above a few GHz, where Faraday rotation can be safely ig¬ 
nored. 

- The new polarization map shows a number of large-scale po¬ 
larization structures, including spurs and filaments. These 
new features are detected clearly only in polarization, in¬ 
dicating a high polarization fraction similar to the brighter 
well-known spurs. 

- The halo of synchrotron emission in the inner Galaxy clearly 
visible at 408 MHz appears to have a low polarization frac¬ 
tion (< 10 %) compared to the high-latitude spurs with typ¬ 
ical polarization fractions of 30-50%. The diffuse high- 
latitude emission is also weakly (though detectably) polar¬ 
ized (<15%). 

- We discuss in detail the characteristics of Loop I. The new 
polarization maps enables us to follow it further into the 
southern Galactic hemisphere, where there is a substantial 
deviation from the small circle fitted to the brighter north¬ 
ern part, presumably due to the inhomogenous environment. 
Radio depolarization. X-ray absorption, and tomographic 
mapping of the local ISM all suggest that the loop is at least 
several times more distant than the traditional 120 pc. 

- The projected magnetic held outside the North Polar spur 
parallels the loop, supporting the model that the loop is in 
part shaped by its interaction with the ambient magnetic 
held. 
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- The filaments projected inside Loop I are mostly associated 
with it, possibly a system of internal shock waves created 
when a re-energizing supernova shock reflected from the old 
boundary of the cavity. Alternatively, in accordance with the 
model of Vidal et al. (2015), they might be a system of “il¬ 
luminated” field lines on the back hemisphere of the loop 
cavity. 

“ On several grounds we are not convinced by the suggestion 
of Liu et al. (2014) that emission from the loop is seriously 
contaminating derived CMB maps. 

- We trace the magnetic field in Loop II for the first time. It fol¬ 
lows the path of the loop, as in the other examples. We show 
that the most prominent feature in the original definition of 
Loop II is actually unrelated, and so the loop is somewhat 
smaller and fainter than previously believed. 

- The proposed extension to Loop III in the south suggested by 
Vidal et al. (2015) is largely an artefact in the WMAP maps. 

- The South Polar Spur, previously identified as the low-/ edge 
of Loop II, is actually the high-/ edge of a smaller struc¬ 
ture. Just like for Loop I, there is a filament of cold material 
(traced by H i and thermal dust) running just outside it, with 
parallel magnetic field. From the H i velocity the distance is 
at most a few hundred pc. 

- Loops V and VI (Milogradov-Turin & Urosevic 1997) and 
SI (Wolleben 2007) are not physical structures; our map 
shows that major components of the proposed loops devi¬ 
ate radically from the expected path. These components are 
independent features, with much smaller angular size than 
the suggested loops. 

- We find two examples of smaller and more distant loops, one 
of which arcs around the Cyg X star-forming region. Both 
are associated with polarized dust features. 

- We observe a clear outline of narrow highly polarized in¬ 
tensity around the northern Fermi bubble above the Galactic 
centre. The spectral index appears to be much fiatter {p = 
-2.54 + 0.16) than other filaments, but similar to that found 
for the microwave haze. This is indicative of a more ener¬ 
getic population of electrons along these lines of sight. We 
note that the Southern extension of the Fermi bubble extends 
further south than Loop I, which probably rules out an asso¬ 
ciation of Loop I with the Fermi bubbles. 

- We have discovered examples of anti-correlation of Ha and 
polarized synchrotron intensity. The most striking is a high 
latitude Ha filament that appears as a local minimum in 
polarization. Because of the low density of ionized gase 
we can rule out Faraday depolarization. The observed anti¬ 
correlation appears to be real in that there is no synchrotron 
emission emanating from the region of the Ha filament. The 
negative velocity of Ha and the longitude imply that the fil¬ 
ament could be lying at a distance of at least 2 kpc. 

- We correlated the polarized maps with the Commander AME 
map. We found that there was no significant correlation 
except for the Pegasus plume region, where a low level 
(9.0 ± 1.9 %) of polarization was observed, consistent with 
contamination from surrounding synchrotron radiation. 

Although we are beginning to understand the diffuse emis¬ 
sion at Planck frequencies, there is still a long way to go. 

In intensity, the component separation is still a major difficulty. 
In particular, we have not been able to constrain the synchrotron 
spectral dependence, due to the degeneracy and confusion with 
free-free, AME, and CMB. New data in the frequency range 
2-15 GHz are needed to more fully characterize all of the low- 
frequency foreground components. Experiments such as S-PASS 


at 2.3 GHz (Carretti 2011), C-BASS at 5 GHz (King et al. 2014), 
and QUIJOTE at 11-19 GHz (Rubino-Martm et al. 2012b) will 
improve the situation considerably. 

In polarization, the emerging picture is that synchrotron 
emission is by far the dominant polarized component below 
50 GHz, with very little contribution from free-free and AME. 
The diffuse synchrotron emission has now been measured with 
good S/N ratio for a large fraction of sky. However, there are 
still regions where the S/N ratio is low. More importantly, the 
detailed spectral dependence of synchrotron emission is still un¬ 
known. 
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Appendix A: Simple ILC method 


dust. We omit maps with frequencies in the range 50-120 GHz, 
where the foreground components are weak and their spectrum 
(especially the dust components) is not well defined. The weak¬ 
est point of our modelling is the assumption of a fixed spectral 
shape for the thermal dust, but this component is very weak be¬ 
low 50 GHz so that errors in the dust spectrum there have little 
effect; the spectrum is mainly used in the interpolation from 353 
to 143 GHz needed to generate a dust-cleaned CMB template. 

We solve Eqs. A.2 and A.3 by restricting the analysis to n 
frequency maps, where n components are dealt with, to make 
Mki square, and then 


w = M 


( 0 


4^03 


(A.4) 


In our final application, we have more bands than constraints 
and so we found the weights by maximizing the S/N ratio for our 
reference component: 

= , (A.5) 


The well-known ILC methods return an image of a component 
with a known spectrum in the presence of other components 
whose spectra are either not explicitly considered (Bennett et al. 
2003) or only known in a subset of cases (Remazeilles et al. 
201 la,b). In Sect. 3 we have the more limited aim of eliminating 
one or two of the major components with known spectra at low 
frequencies (CMB and free-free emission), and returning a map 
reasonably close to the sky at 28.4 GHz with these components 
omitted. 

We model the maps as sums of components with spatially- 
invariant and known spectra, plus noise and other unmodelled 
foreground emission: 

tjk = MkiAij -r Njk + Ojk (A.l) 

where the indices are i for component, j for pixel and k for band. 
Aij is the amplitude. The mixing matrix Mki = fi(yk)/UkCik, 
where ftivk) is the spectral form in Rayleigh-Jeans brightness 
temperature at the band reference frequency Vk, Uk is the unit 
conversion from Rayleigh-Jeans to thermodynamic temperature, 
and Cik is the colour correction (Planck Collaboration IX 2014; 
Planck Collaboration II 2016).^^ For the CMB, component 0, we 
have Mko = 1. 

Our aim is to find the weights for the frequency maps that 
satisfy the conditions 

J^MkiWk = 0 (A.2) 

k 

for the components to be eliminated, subject to the normalization 
condition 

^ Alk?>Wk = Mo3 (A.3) 

k 

where band ^ = 0 is 30 GHz and component i = 3 corresponds to 
3. p = -3 power law, roughly appropriate for synchrotron emis¬ 
sion, but also quite close to the effective spectrum of spinning 

Colour corrections and unit conversions for WMAP were derived 
in the same way as for the LFI, using the released WMAP bandpasses. 
Colour corrections were based on the spectral index of the model eval¬ 
uated at Vk, except that we used the explicit dust colour corrections pro¬ 
vided by the UC_CC code for the HFI bands (Planck Collaboration IX 
2014). 


where, to keep the weights spatially uniform, we take cr^ as the 
median pixel variance for band k. We split the bands into a set k' 
whose weights are varied by the fitting routine (AMOEBA in IDL), 
and a set k whose weights are found by linear regression. Our 
constraint is now written 

^ Mki^^k = ^ Alk'iw[. (A.6) 

k k' 

Each AMOEBA iteration provides a trial set of values that are 
used to find Wk via 


Wk 


^ ^kl ^i3^03 - ^ Mk'iWk' 
/V k' > 


(A.7) 


and the full set of weights is then used to evaluate the S/N ratio. 
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